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SECTION 1 


INTRODUCTION and SUMMARY 

This report is a description of the development effort for the MAESTRO program 
to be used for the in-flight decision-making process during the tr an s lunar and 
lunar orbit adjust phases of the flight of the Radio Astronomy Explorer - B 
(RAE-B). MAESTRO is an acronym for Mission Analysis Evaluation and Space 
Trajectory Operations. The program serves two functions, performance and 
evaluation of preflight mission analysis and in-flight support for the midcourse 
and l una r orbit insertion command decisions that must be made by the flight 
director. 

The all- FORTRAN program was designed to be a flexible mission analysis tool 
that can be readily modified to analyze and support different types of spaceflight 
missions. One of the principal guidelines in the MAESTRO development effort 
has been to reach an equitable balance between efficiency and flexibility. The 
advantage of this type of development was demonstrated in the summer of 1972 
when, within a few months, the program was modified to support the IMP-H 
mission analysis and in-flight operations. 

This document supplements a previous report on the mission control system 
design study ( Final Report for NAS5-11796). The skeleton system developed 
in the earlier contract is described there at length along with the analyses that 
helped set the criteria for development of the actual working program. 

The development of the system has been divided into three basic tasks as follows: 

I System Integration 

n On-line Control and Displays 

III Operational Logistics and Testing. 

The report is arranged according to these tasks except that the extensive midcourse 
guidance analysis has been included as a separate section. 
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Section 2 contains a basic description of the program capabilities along with 
discussions of the analyses that accompanied the development of those capabilities. 
The program options are divided into operational modes that are invoked by the 
setting of a single control flag. Each of the options uses the basic trajectory 
propagation software that forms the heart of the program. This trajectory 
propagation software itself contains several options for selection of the trajectory 
propagation method, numerical integration technique or analytic propagation 
technique. The program options include straightforward trajectory propagation or 
ephemeris generation, orbit insertion retro- or kick-motor firing analysis, mid- 
course guidance analysis with any of a number of guidance laws and targeting 
constraints and a fully integrated midcourse verification mode. The program 
contains orbit prediction capability in two forms of different speed and accuracy, 
a post injection .trim analysis with several levels of sophistication, along with 
shadowing, occultation and tracking visibility calculations as options. These 
capabilities are not mutually exclusive and the options can be run in tandem to 
form the Monte Carlo mode that provides statistical information on expected fuel 
requirements, orbit trim requirements, spin-axis-sun angles at maneuvers, 
distributions of execution errors and a great deal of more detailed information 
related to the specific goals of the RAE-B mission. 

The Monte Carlo capability can be run either as a production type, pre-launch 
mission analysis tool or on-line during the flight to help the flight director make 
decisions on the basis of maximizing the probability of mission success. These 
capabilities and the available optional features are described in detail in 
Section 2. 

Section 3 contains an extensive discussion of the midcourse guidance capabilities 
of MAESTRO. Detailed descriptions, with explicit examples, are given for the 
various quidance laws. These laws, and the entire guidance philosophy, have 
evolved through a combination of theoretical and numerical experiments aimed 
at providing flexibility and speed of computation. 

The midcourse quidance calculations are normally accomplished with the 
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Multiconic trajectory propagation technique because of its speed and flexibility. 

These calculations can, however, be performed with full numerical integration 
if desired. The guidance end-point constraint variations are evaluated by the 
secant method and the required mid, course velocity is obtained by the usual 
linear analysis in each iteration of a numerical search procedure. 

The midcourse quidanee analysis features an analytic pre-targeting procedure 
that is based upon a modified sphere of influence patched conic technique and 
internal gradient evaluation control logic that greatly reduces the required 
iterations in the midcourse analysia 

The results of the quidanee calculations are further used in the evaluation of 
expected second midcourse fuel requirements, expected errors in the lunar approach 
trajectory, and the expected variations in the post-injection lunar orbit. These 
quantities are available, on option, to the flight director to aid in the selection of 
the best overall midcourse maneuver time and attitude. The guidance analysis 
mode provides for selection of minimum fuel, fixed time of arrival, fixed 
pre-retro selenocentric energy, fixed post-retro energy, variable target 
inclination procedures and a variety of overburn or underburn options for 
these guidance modes. 

Section 4 is a description of the methods for on-line control and the printed 
displays of the MAESTRO program. The program includes an input editor 
that assigns preset values to the control variables according to the setting of 
the mode flag. These preset values can be changed with little effort but the 
objective of the automatic input scheme is to permit the program operation to 
be semi-automatic. The preset inputs have been arranged in a hierarchy so 
that certain analyses can be run by the setting of the single control variable. 

This scheme greatly reduces the chances of human error in operation, especially 
under conditions of stress during an actual flight. 

Also in Section 4 are brief descriptions of the printed displays that are used to 
form a director's book for reference and presentation of results. These displays 
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are designed to be comprehensive but also readily understandable. Many of the 
displays are split apart with the bulk of the information printed separately from 
a concise summary sheet that presents a readily- grasped evaluation of the situation. 
Several examples of the displays are included in this report and the reader is 
referred to separate documents for additional examples and explanation* 

Section 5 contains explanations of the methods of using the program for in-flight 
operations support. The first subsection includes a typical flight profile 
outline along with an approximate time-line for the in-flight control operations. 

The hypothetical flight pro file is only one of several similar possibilities 

but the outline contains all the operations that might be required in a "worst case" 

situation. The discussion emphasizes the importance of a close man-computer 

relationship during the in-flight operations and the need for a smooth flow of 

information between the separate parts of the overall control system is made 

obvious. 

A second subsection contains discussions of possible failure modes with fairly 
detailed outlines of three of what appear to be the most likely contingineies. 

The areas of possible failure are catagorized in such a way that the reader 
can associate his own specialty with one of the modes. The presentation is 
intended to stimulate the reader's imagination so that possibilities that may not 
have been considered will be brought to light. Failures are classified broadly as 
1) Hardware malfunctions, 2) Computer Software Errors, and 3) Human errors. 

In addition to this generalized classification of possible failures, the three 
specific areas of most import are suggested to be: 1) Midcourse and trim 
propulsion sytem failure, 2) Attitude control and determination system 
malfunction, and 3) Retro rocket failure. These three possibilities are 
discussed in some detail. 

A third subsection contains descriptions of the system testing procedures, the 
simulations that have been held to date, and recommendations for the extensive 
testing that is to take place during the next few months. Several of our early 
mistakes are discussed and suggestions are made for minimizing the possibility 
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of such mistakes in the future. Emphasis is placed upon the need for complete 
and rapid communication between the components of the mission control system 
and upon the need for a close man-computer relationship during the in-flight 
operations. Also stressed is the need for a smooth flow of information between 
the three major components of the computer software system. 

A closing section includes a brief summary of the important aspects of the 
MAESTRO system development and operation. A short list of recommendations—. 
is presented along with some suggestions for the future development, use, and 
upkeep of the MAESTRO program. The reader is reminded that this document 
is complemented by four other reports: the MAESTRO Programmer's Manual 
Reference (6), the design study final report Reference (1), the MAESTRO User's 
Manual Reference (2), and the Parameter Estimation Notebook Reference (7) 
which might be considered an appendix to this report. The non-technical reader 
who wishes to skip to the conclusions and recommendations should first examine 
the operations flow diagram of figure 4.1. 



SECTION 2 


SYSTEM DEVELOPMENT 

The MAESTRO system is an outgrowth of the skeleton system described in 
Reference 1. The skeleton system demonstrated the feasibility of developing 
a mission control system that has the capability of propagating different kinds 
of trajectories as efficiently as possible. The trajectory propagator in the 
skeleton system was developed using the philosophy that no one trajectory 
propagation technique is most efficient for all trajectories. For example', 
it was found that Cowell was most efficient for a nearly rectilinear orbit and 
Encke’s method is most efficient in other situations. Thus, the trajectory 
propagator in the skeleton system had various trajectory propagation techniques 
and numerical integration techniques. MAESTRO has the same trajectory 
propagation tool as the primary workhorse in the system. There are several 
mission analysis modes built around the trajectory propagator. These modes 
are: 

1. trajectory propagator 

2 . retro motor firing analysis 

3. midcourse analysis 

4. Monte Carlo analysis 

5. midcourse verificiation analysis 

6. orbit evolution 

7. post injection trim analysis 

The following sections describe each of these modes. 

It is important to note that MAESTRO was not designed to be restricted to a 
particular mission. Although many of the displays are concerned with the mission 
objectives of the RAE-B mission, the entire system can be easily adapted to inter- 
planetary or Earth orbit missions. In the Summer of 1972, MAESTRO was quickly 
modified and implemented as the mission control program for the IMP-H mission. 
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2 . 1 TRAJECTORY PROPAGATION MODE 


In this mode MAESTRO is used as a three-degree-of- freedom trajectory 
simulation program with the capability of propagating the trajectory with any 
of the following trajectory simulation schemes: 

a. Cowell's method 

b. Encke's method 

c. numerical integration of the classical elements 
with true anomaly as the fast variable 

d. numerical integration of the classical elements 
with mean anomaly as the fast variable 

e. averaging of equations in T, d" 

f. multiconic 

g. numerical integration modified elements 

h. averaging of equations in "g" 

The first four trajectory propagation schemes are commonly used techniques, 

A description of how they were incorporated into MAESTRO is described in 
Reference 1. The averaging used in schemes "e ,T and "h" are detailed in 
Section 2 . 6 of this report. Multiconic is an approximate method which gives 
very good results for Earth-Moon trajectories. Its chief advantage is that 
it is extremely fast. A description of this technique was also presented in 
Reference 1. Scheme "g" is similar to the integration of the classical elements 
except that the elements used do not have a discontinunity as the eccentricity 
approaches zero. The terms to be integrated are: 

1. semilatus rectum, p 

t 

2 . e sin cj 
V 

3 . e cos <o 

4. &> + f 

5. inclination, i 

6. longitude of the ascending node, 

where e is the eccentricity 

to is the argument of perigee 

f is the true anomaly 
7 



The derivatives of the above quantities are numerically integrated to obtain 
the instantaneous orbital elements. The equations which define these 
derivatives are presented in Appendix A. 

The user also has the option of selecting any one of a variety of numerical 
integration schemes. The schemes available are 

a. Ten-cycle seventh-order Runge-Kutta 

b. Twelfth-order multistep 

These schemes were described in Reference 1. Reference 1 also presented 
an analysis to determine the relative accuracy of these schemes. In that 
report it was stated that the seventh-order Runge-Kutta appeared "best.” 

This integration scheme is still considered the best with most of the trajectory 
propagators; however, it was found that Cowell can achieve greater accuracy 
in the same computing time when used with the twelfth-order multistep. The 
twelfth multistep method possess the limitation that it cannot be easily used 
with a varying compute interval, and requires a cumbersome startup procedure. 
The other capabilities of the trajectory propagator in MAESTRO are 

1. Planetary Epheme rides 

a. DE 69 tape read for all planets 

b. mean elements for Sun 

c. mean elements for Moon 

d. o sculating elements for Moon 

e. DE69 disk read utilizing direct read feature on IBM computer 

2. Lunar Gravity Models 

a. Ll field 

b. JPL - 15 by 8 field 

(a high-order geopotential can be included in the 
same logic for close-Earth oi'bit analysis) 

c. Post Mariner 9 Mars field 

3. Input Coordinate Systems 

a. mean equator and equinox of 1950 

b. mean equator and equinox of date 
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c. mean equator and ecliptic of date 

d. true equator and prime meridian 

e. true equator and equinox of date 

f. true equator and node 

g„ spherical Earth centered 

4. Output Coordinate Systems 

a „ mean equinox and ecliptic of date 

b. true equator and prime meridian 

c. mean equator and equinox of 1950 

d. true equator and equinox of date 

5. Rocket motor firing with variable thrust/weight flow capability 

6. Transfer orbit and in-orbit shadow calculations for both umbral 
and penumbral passages 

7. Solar pressure 

8. Run termination at 

a. input final time 

b. closest approach to a target planet 

c. farthest distance from a target planet 

9. Compute interval 

a. preset table with up to ten values as a function of time 
b„ automatic 

10„ Variable printing logic 
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2 . 2 RETRO MOTOR FIRING ANALYSIS 

The retro motor firing analysis or approach analysis is used to determine 
the attitude and time to fire a retro motor to achieve desired mission constraints. 
In this analysis it is assumed that all midcourse maneuvers have already been 
completed, thus the approach trajectory (trajectory before retro firing) is 
fixed. The analysis can be performed in either of two levels of accuracy. 

In the most accurate level, all segments of the approach trajectory are 
numerically integrated and retro motor is also numerically integrated. In 
the lower accuracy level the orbital elements of the approach trajectory are 
determined. The state at each firing time is determined by stepping along this 
constant set of elements. The retro motor is not numerically integrated; 
instead the impulsive velocity of the motor is used. The orbital elements of 
the approach trajectory are determined by numerically integrating the approach 
trajectory to the desired stop time (or closest approach if the appropriate 
inputs are set). The orbital elements at this time are used as the constant 
elements of the approach trajectory. 

No matter which level of accuracy is used, the analysis is performed in the 
same manner. Firings of the retro motor are simulated along the approach 
trajectory at a series of attitudes. The mission constraints are determined 
for each firing at each attitude and displays are presented so that the user can 
rapidly determine the firing time and attitude which best satisfy the mission 
constraints. The range of firing times and attitudes are determined in the 
following manner: 

1. Firings are made between input range of true anomaly about 
the stop time of the approach trajectory. The attitude range 

is an input cone angle about the velocity vector at the stop time. 

2. Firings are made from asymptote to asymptote on the approach 
hyperbola. The attitude range is input as initial right ascension 
and declination, increment in right ascension and declination, 
and the number of right ascensions and declinations to be tried. 

The approach trajectory must be hyperbolic. 
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3. Firings are made between input times. The attitude 
range is input as in 2. 

Although the displays currently available are designed for the RAE-B mission, 
displays could easily be changed to accommodate other missions. In this 
mission the most important mission criteria is achieving the desired Lunar 
orbit with sufficient maneuver fuel. Thus, a display is presented to rapidly 
access the retro firing attitude which l’esults in the minimum maneuver 
fuel. Other displays are presented which have more complete information 
about the firings made at each of the attitudes. Samples of these displays 
are presented in the useds guide, reference 2. 
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2.3 


MIDCOURSE ANALYSES 


The midcourse analysis capability in MAESTRO is described in Section 3. 
This section contains the theory, run setup and practical results. 
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2.4 MONTE CARLO ANALYSIS 

The purpose of the Monte Carlo analysis is to determine the probability of 
mission success. The probability of mission success is determined by 
"flying” the mission numerous times with errors applied to the initial state 
and to each of the maneuvers. Each of the samples are analyzed to determine 
whether they satisfy mission constraints. The probability of satisfying a 
particular constraint is determined from the sum of the number of samples 
which satisfy the constraint divided by the total number of samples tried. 

A typical mission flight plan could consist of the following: 

1. Apply tracking errors to the initial state and fly to the 
first midcourse execution time. 

2. Determine the first midcourse correction, apply errors to 
the correction and fly to the second midcourse execution time. 

3. Determine the second midcourse correction, apply errors to 
the correction and fly to the Moon. 

4. Determine the retro motors firing time and attitude. Apply 
errors to the retro and burn retro. 

5. Determine if final Lunar orbit satisfies mission constraints. 

The analysis currently has the capability of simulating a flight plan similar 
to the one outlined above or a plan with only one midcourse and retro or only 
a retro maneuver. The Monte Carlo analysis also has the capability to determine 
the first maneuver from the current state. This maneuver, with errors, is 
applied to all of the samples. This capability should be used when the errors 
in the estimation of the initial state approach a constant value. When this point 
is reached the state will not be known to a more accurate level in the future. 
Thus, the first maneuver must be determined using current state and not the 
current state with estimation errors. 

The model used to simulate errors in the initial state is a covariance matrix. 
This matrix is calculated by the oribit determination program and passed to 
MAESTRO, or can be input. 
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itrror Models 


The tracking error covariance matrix, P, is a 6 x 6 positive definite matrix of 
anchor vector estimation errors. 

P = E (x - x) (x - x)* (x = anchor vector, x = true state) 

There is a preferred coordinate system in which components of the error factor, 
y, are uncorrelated. The (similarity) transformation between y and x - x is S, 
an orthogonal matrix which diagonalizes P. 

x - x = Sy 

The diagonal matrix, D, defined by 
0 = 8* PS, 

has as its diagonal elements the variances of the uncorrelated components of y. 
Scaling a white noise (uncorrelated random numbers) sampled by the standard 
deviations of y, we obtain an error vector in uncorrelated coordinates with com- 
ponents, y^. 

y. = d i n (0, 1) 

(dj is the i-th diagonal element of D, n (0, 1) is a sample random number from a 
normal (Gaussian) distribution of mean 0 and variance 1.) A sampled errant 
trajectory state is then computed from 

x = x + Sy. 

Execution errors for the midcourse and retro maneuvers are computed from 
assigned standard deviations of pointing and velocity errors. The velocity 
error is treated as a proportional error normally distributed plus a resolution 
error uniformly distributed. 

<s v = vtfyn (0,1) + (. 0001 km/sec) u 
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(v is the velocity impulse magnitude and u is a random number from a distri- 
bution uniform on the interval v{ 5<u < .5 | ). The retro velocity error is 
formulated without a resolution error. The pointing error is formulated as 
two independent errors normally distributed along mutually orthogonal axes 
which are both orthogonal to the maneuver impulse direction. Let the maneuver 
impulse be denoted by V and its direction by V. If K is the unit polar axis of 
V's coordinate reference frame, we can construct unit vectors normal to V. 



N = V x E 


K - V (K« V) 
cos (dec) 


The "eastward" and "northward" pointing errors, Q e and , are computed by 
scaling random numbers with the input pointing errors, • They are then 
converted to velocity deviations using the small angle approximation. 


f e = vo a n x (0, 1) 


f n = vo a n 2 (0,l) 

The misdirected velocity vector has the direction 
V' = (vV + <e E +c n )N /Vv 2 + ( q Lr a 2 


And the magnitude v =« v , so that it may be programmed as follows, 


V' = v +6 v 


V 2 2 2 

1 +<r a (nj + n 2 ) 
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The Monte Carlo analysis demonstrated the necessity to develop some method 
to automatically determine the retro motor's firing time and attitude. In 
earlier analysis models, the retro was simulated at closest approach to the 
Moon opposite to the velocity vector. This firing strategy was sufficient when 
the spacecraft passed close to the aiming point at the Moon. However, small 
midcourse execution errors will cause relatively large errors at the Moon, 

Thus, some other retro firing strategy had to be determined. The approach 
finally selected used a steepest ascent optimization procedure to determine the 
attitude which minimizes a particular function. The function minimized is the 
minimum post injection trim velocity on a Lunar approach hyperbola. This 
velocity is determined by firing the retro at various points along the approach 
hyperbola and calculating the trim velocity to achieve the desired final Lunar 
orbit. The minimum of all of the calculated trim velocities is the value used 
in the function to be minimized. 

The trim philosophy is to perform a two - impulse Hohmann transfer to achieve the 
desired circular orbit and to make a third plane -change maneuver at the 
largest node of the intermediate orbit. Whether or not this new policy is 
optimal is difficult to decide. It is certainly not "fuel optimal," as the litera- 
ture shows, for large plane change situations. In the normal case of small 
plane change, the fuel optimal policy will probably yield fuel savings that are 
inconsequential in comparison with the complexity and risk of performing 
additional attitude reorientations and engine ignitions. In any case, the logic 
for determining the required trim fuel, given a specific post-retro orbit, is 
an isolated subroutine which can be easily modified as more sophisticated 
trim maneuver logic is incorporated. The point important to the above dis- 
cussion is that a value can be determined that represents the minimum fuel 
required to change the orbit from its post-retro size and orientation to the 
desired lunar orbit suitable for long-term antenna stability. 

A more detailed description of the steepest ascent optimization procedure is 
shown in Appendix B. 
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The Monte Carlo output presents the running means and standard deviations 
of the variables of interest. These variables include the fuel requirement and 
spin axis - sun angle for each maneuver as well as the post- retro orbital 
elements. A statistical summary is printed at the end of the run. This 
summary contains the mean, standard deviation, maximum value encountered 
and minimum value encountered of the quantities mentioned above. The analysis 
also has the provision to output starting conditions for each sample so that any 
particular sample can be studied in detail. 



2.5 


MIDCOURSE VERIFICATION ANALYSIS 


The midcourse verification, analysis was included in MAESTRO in order to 
provide a check of other analysis and provide data useful for real time decision 
making. The mode can be used to verify any midcourse motor firing. This 
includes firing for midcourse corrections on the transfer trajectory or firings 
in orbit for post injection trim. 

In this analysis, the percise numerical integration techniques are used to 
propagate the state from the initial time to motor ignition. Next, the motor 
is numerically simulated and the spacecraft propagated to closest approach 
to the target planet if desired. The final time can be set to the engine burn- 
out time so that the analysis ends there. This would probably be desired for 
post injection trim verification. 

A doppler analysis of the motor burn is also presented. The output from 
this analysis provides histories of the nominal spacecraft thrust, weight 
and the velocity away from the visible tracking stations assuming nominal thrust 
characteristics. The velocity away from the tracking stations can be used 
to calculate the doppler shift expected during the motor firing. The expected 
doppler shift can be compared to the actual shift experienced in real-time and 
correction actions can be initiated if a non-nominal burn is indicated. 

The midcourse verification analysis also presents a list of tracking station 
elevation and azimuth angles at motor ignition time. Thus, tracking station 
coverage is automatically available for any motor ignition time by the use of 
this mode. ‘ 

The post maneuver orbit is presented at closest approach to the target planet 
(or final time). This orbit should be compared to the orbit predicted by the 
midcourse analysis or post injection trim analysis. The orbits should be 
the same if same techniques were simulated in both analyses. For example, 
if impulsive velocity was assumed in the midcourse analysis and not used 
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in the verification run, differences arrising from the impulsive velocity 
approximation will appear when comparing the two runs. 

The displays obtained from this mode are shown in the user's guide. 



2.6 


LUNAR ORBIT PREDICTION - Numerical A . ve raging 


One of the primary objectives of the RAE-B mission is to establish a lunar orbit 
whose eccentricity will remain near zero for at least a year. The complexity of 
the lunar gravity field makes it difficult to perform approximate predictions of 
the orbit evolution by the usual analytic techniques. Because of this difficulty 
and because of the need to predict the long-term motion of the orbit, a numerical 
averaging technique has been developed and implemented in the MAESTRO program. 
This technique was described briefly in Reference (1) and more thoroughly in 
Reference (5). The following subsections contain discussions of the averaging 
method, a newly implemented numerical startup procedure to obtain initial mean 
elements, and some examples of the method's efficiency in the RAE-B orbit 
prediction problem, 

2.6.1 The Averaging Technique 

The method as implemented in MAESTRO is described in detail in Reference ( 5 ). 
The essence of the technique is the averaging of the planetary equations (see 
Appendix A) by numerical quadrature over one revolution of the spacecraft in its 
orbit. It is shown in Reference ( 5 ) that the system of ordinary differential 
equations describing the osculating orbital elements can be transformed, by the 
averaging operation, into a system of ordinary differential equations describing 
the mean values of the orbital elements. 

Let the system of planetary equations given in Appendix A be represented by 

E. = F. (E., t) i=l,6; 3=1,6. 

x 1 j 

The time average of this system is defined as 
_ t+r/2 

E=— / F. (E., t)dt i=l,6;j=l,6, 

T I 1 J 

t-r/2 

where r is the orbital period. 
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The integration variable is transformed to the Keplerian true anomaly, f, as 



2 

n P 

2 IT jU P 


/, 


f(t+r/2) 

r 

f(t-r/2) 


(1 = e cos f ) 


It should be noted that . 

t +r/2 
o 

-~T“ = “ / E (t) dt = E (t +r/2) - E (t -r/2) 

dt T J o 2 

t o - T /2 f. 


and that in terms of similarly defined mean values of the elements. 


d 

dt 


(E) 



1 

T 


t + r/2 
/ E (t) dt 
4-r/2 


dE 

dt 


provided that we maintain the understanding that the orbital elements remain 
constant during the averaging integration. Thus, by these few assumptions we 
have traded the system (1) of ordinary differential equations in terms of the 
osculating elements for the system (4) of ordinary differential equations in 
terms of the mean orbital elements. 


The exact form of the averaging as implemented in MAESTRO is described in 
the MAESTRO programmer’s manual under subroutine AVEQNS and the averaged 
equations are given in Reference ( 5 ). 

2.6.2 Numerical Startup 


For some time there have been discussions of how to start the numerical 
averaging technique (Methods 5 and 8). Past studies have been conducted with 
the assumption that the averaging integrations were started with known values 
of the mean orbital elements. These initial mean elements were assumed to 
be available by virtue of some process whereby they can be related to osculating 
elements at some epoch. This assumption has not been unreasonable because 
of theoretical arguments that every almo st-periodic function has a mean motion. 
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We can assume, then, that the time history of the mean orbital elements is 


approximately the same as the history of the running mean of the osculating 
elements. In practice, we can proceed as if for every set of osculating elements 
at some epoch, there is a corresponding set of mean elements representing the 
average values of the osculating elements over the almo st-period of the osculating 
motion. The problem, in the real world, is to relate the osculating elements 
to the mean elements at some epoch. 


In Reference 5, it is suggested that the relationship between mean and osculating 
elements can be obtained by performing a one -revolution integration of the 
osculating orbital elements and by simultaneously forming the integrals 
t 

E. = —— f E. dt, i *1,6, 

i t J l 

0 

where E. represents the osculating elements. Now when the Keplerian period given 
by the running mean semi-major axis is equal to the current time, t, we can obtain 
a mathematically and dynamically consistent one-to-one relationship between the 
mean elements at t/2 and the osculating elements at t = 0. 


This suggestion has been incorporated into the MAESTRO program and the new 
option can be invoked by setting location 1098 to 1. The new option allows one to 
input osculating elements at some epoch and to obt ain a numerically averaged 
solution for the subsequent motion that is dynamically consistent with the initial 
osculating elements. 

Figure 2. 1 shows the time history of the selenocentric argument of pericynthion and 
includes a comparison of the osculating and averaged values of that quantity. The 
motion of the mean elements through the center of the envelopes of the osculating 
elements is very encouraging. The figure is representative of all the orbital 
elements and success of the startup procedure implies a better representation 
of the long-term motion in the realistic case where we must use osculating elements 
to start the numerical averaging technique. 
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NUMERICAL AVERAGING WITH AUTOMATIC STARTUP 


Figure 2 . 1 




2. 6.3 Comparison of Averaging with Actual Solutions 

In Reference (5 ), several examples are given that indicate the speed, accuracy, 
and flexibility of the numerical averaging. It is carefully pointed out that the 
averaging is an intermediate technique that lies, in speed and accuracy, somewhere 
between full numerical integration and the usual analytic techniques. The principal 
advantage of the method is that all kinds of forces can be included in the dynamic 
model without the need for analytic approximation of the disturbing forces. 

Throughout the MAESTRO development effort, experiments have been conducted 
to determine the ’'best" combination of intervals and ordinates to be used in the 
quadrature integration. For the RAE-B orbit, it has been found that excellent 
results are obtained if the averaging interval (the orbit period) is divided into 
three intervals each employing a 6 point Gaussian quadrature. Figure 2.2 
shows the effect of changing the intervals and ordinates for a typical example of 
the tests conducted. 

Finally, Figure 2.3 gives a 200 day time history of the eccentricity for a 
typical RAE-B orbit. This run was started with mean elements obtained by the 
numerical startup procedure and the agreement of the mean and osculating time 
histories represents an obvious improvement over the example shown in Reference (5) 
where the differences between initial osculating and mean elements were ignored. 

The accuracy and timing shown on the figure are typical of the performance of the 
averaging for RAE-B type orbits. 
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Figure 2.3 



POST INJECTION TRIM ANALYSIS 


2.3 


The post injection trim analysis is required to determine the maneuver(s) 
necessary to achieve a desired orbit. The orbit after retro will probably 
not be the desired orbit because the fixed retro velocity could not yield the 
desired orbit, or, if it could, the errors in the retro maneuver itself could 
yield substantial errors. The desired orbit is defined as a circular orbit with 
a specified radius and inclination. The position in the orbit or the argument 
of the ascending n ode is not important. 

There has been many attempts to solve this complex problem analytically. 
Thus, a literature search was initiated to determine if any previous work could 
be adapted to our problem. The search was conducted by Dr. Mary Payne, 
reference 3. Her work indicated that there were a few techniques available 
that may solve our problem. Some of the techniques she suggested were 
techniques that utilize many impulses. It was decided not to use these 
techniques because the many motor firings and attitude changes would degrade 
the reliability without any large gain in system performance* 

The work by T. Sun, reference 4 , seemed the most promising. He developed 
the optimum solution to the two impulse 1800 transfer between non-coplanon 
orbits. It is not necessarily true that the 180° transfer is optimum, however, 
it is optimum in certain planar cases (Hohmann transfer). Since the inclina- 
tion changes in the orbit must be small because of system capabilities, it 
was concluded that the 180° transfer is near the optimum solution. Also, a 
two impulse maneuver is desirable from a reliability standpoint. The follow- 
ing paragraphs present a synopsis of Sun's method and describes how it was 
implemented. 
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Sun’s method determines the optimal two impulse 180° transfer between non- 

o 

coplanar orbits. Since a 180" transfer is specified, the first impulse must be 
applied at the intersection of the initial and final orbit planes. Thus, the angle 
between the initial and final orbit planes and the position on the initial orbit 
where the maneuver is made can be obtained from the spherical trigoinetric 
relationships, see figure below: 


Z 



The angle from the reference plane to the common line of nodes in the initial 
orbit, \ , can be determined from the input initial true anamoly, f, and the 
argument of the ascending node of the initial orbit as, 

X = f ^ w (1) 
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Then, the angle between the two orbit planes, o , can be determined from. 


sin /3 = sin ij sin i / sin \ 
tan -2 = sin ) 


( 2 ) 


sin tan • 


) 


where i is the inclination 


£ is the angle from the reference plane to the common line of 
nodes in the final orbit plane. 

and the subscripts i and d refer* to the initial and desired orbits, 

respectively. 


The radius and velocity components can be determined from standard orbital 
relationships at the initial true anomaly, f. 


The orientation of the transfer plane with respect to the initial and final orbit 
planes is described in the figure below:, 



The angles and nt, describe the orientation of the transfer plane with respect to 
the initial and final orbit planes, respectively. 


If the inclination of the transfer plane with respect to the initial plane is specified, 
then the optimal velocity can now be determined using Sun's equation 10. 


AV = 
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where ^ is the gravitational constant 
r is the radius at f = X 

V rl , V j are the radial and transversal velocity components on the 

initial orbit relative to local velocity ( V=^/ jU /r ) 

n is r ,/ r. 
d 1 

The equation above is somewhat simplified from Sun's equation since the final orbit 
is circular. Thus, 

V r2 ’ 0 
V T2 = + 

If the velocity of each trim maneuver is desired, then 
AY 

- — 1 H 

Av V 5 ” 

AV„ 


(4) 
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cos 
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(5) 


where AV^ and A V^, denote the magnitudes of the first and second trim 
maneuvers. 

The direction these impulses are applied can be determined by noting the following 
relationships. 

AV “V -V 
R1 RT1 R1 

"'' V NI " " V T1 Sln “l 

AV T1 = V TT1 ’ V T1 C ° S <4 
The second trim is determined in a similar manner as 


( 6 ) 


AV „ 

= V 

- V 

r2 

R2 

RT2 

AV N2 

= -V 

T2 

Sin 

4V T2 

= V 

T2 

c °s OJg " 


( 7 ) 


TT2 
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In the above equations the components of velocity are defined as follows, 


V . 

Hl,2 

V 

Tl,2 

V 

TT1.2 

V 

RTl, 2 


radial component of velocity of the 
initial or final orbit. 

Transversal component of velocity of the initial or 
final 1 orbit, in the initial or final orbit plane. 

Transversal component of velocity of the initial or 
final orbit written in the transfer plane. 

Radial component of velocity of the transfer orbit at 
initial and final orbit crossings. 


All components of velocity except the radial components of the transfer orbit in 
equations 6 and 7 are fixed by specifying the initial and final orbits. The radial 
component of velocity is determined from the condition that the total trim velocity 
is to be minimized. The total trim velocity is 


AV = £V 


j+AV 2 = 4 


'* V H 2 +iV N +(V rTl- V rl )2 




+ VaV t2 + AV n2 0^ r2 “ V rT2^ 


Also, 

V = - V 
rTl rT2 


since a 180° transfer is specified and the final orbit is to be circular. Now, the 
partial derivative of the trim velocity with respect to V rTl can be written as, 


( 8 ) 


( 9 ) 


dV 


rTl 


V - v 
■ xXJ Ixl_ 

AV. 


V 

+ . rTl 

AV 0 


= 0 


or 


rTl 


AV* YrL 


( 10 ) 
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The components of the trim velocity obtained from equations 6, 7 and 10 describe 
the trim velocity with respect to the transfer plane. The trim velocity vector in 
the same coordinate system of the initial orbit is obtained through a three Euler 
angle rotation pictured in the figure below. 


Z 



In the figure above, x corresponds to the radial direction, y to the transversal 
direction and z to the normal direction. 

The angular elements of the reference orbit define the Euler angles. Thus the 
transformation from the x, y, z system to the X, Y, Z system is 



cos ib cos O 
-cos i sin Osin \b 

-sinji cos Q 
-cos i sin Ocos ib 

sin i sin 0 - 


X 1 

■ 

cos i|) sin O 
+cos i cos Osin 0 

-sin ib sin 0 
+cos i cos O cos ib 

-sin i cos 0 


y l 

) 

_ sin i cos ib 

sin i sin 0 

cos i 


z ) 


where & = f + to 


The above equation is used to transform the trim velocity components from an orbit 
plane coordinate system to the system of the reference orbit. 


Optimum Inclination of the Transfer Plane 


The condition for the optimal orientation of the transfer plane is expressed by 

2 1/2 


sm to 2 _ 

sin 




1 - 2/^ cos to 2 +P 2 
1 - 2p x cos oJjl 


z V 

) 


32 



1/2 


where = 


= jfm r d / ri 


r + r 
d 1 


§,Wri/r d 


T1 


vl/2 


r + r 
d i 


Equation 12 along with the condition 

= <? + Wj 

yields a set of equations which can be solved for qj or qv> to yield the optimum 
orientation bf the transfer plane. The solutions to equation 12 resulted in a sixth 
order polynomial in sin cj . The equation was solved numerically in order to avoid 
the cumbersome task of solving a sixth order equation, A Newton-Raphson procedure 
was employed to determine the solution to equation 13. Sun, in the refei“ence, states 
that the solution is unique. Thus, the task of finding multiple solutions with the 
Newton-Kaphson method is not required. 

The optimum 180° transfer solution discussed above determines the maneuver to 
go from a fixed point on the initial orbit to the final orbit. We are not constrained 
to any specific point on the initial orbit, except for spin axis- sun angle considerations. 
Thus, it is necessary to determine the position on the initial orbit to make the first 
maneuver that results in the lowest trim requirement. The optimum position on the 
initial orbit is determined by a search of the entire orbit, A half-interval iteration 
is used to find the local minumum and the position used is the minimum of all the local 
minimum s. 

The post injection trim analysis also determines the Hohmann transfer maneuver 
to correct the eccentricity and raduis. This result is useful as a teference and shows 
how much fuel is expended to change the inclination of the initial orbit. 
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SECTION 3 


i 


MIDCOURSE GUIDANCE CAPABILITY 


3 . 1 Theory 

The constraint error vector, is related to the control vector, U, by the 
non-linear functional expression 

* tf(U). 

The midcourse guidance problem can be stated in the form: "Find a particular 
control vector, U , which nulls the constraint error vector, >£/." We shall now 
proceed to define ^ and'U specifically for the RAE-B mission and to describe the 
implemented solution to the midcourse guidance problem. 

3. 1.1 Controls 

The control vector, U, for the RAE-B midcourse guidance capability is 

identified as the midcourse velocity correction impulse vector, AV . The program 

m 

has capabilities for treating U: 

1. strictly as a velocity impulse applied at midcourse ignition time and 

2. as a set of numbers to be transformed to obtain burn duration and thrust 
direction which will render the post-burn state through integration of the 
powered equations of motion. 

The transformation of U for finite burns is 

6w = w ( 1 - e U/ ^ C ) 
o ' 

^ = f(6w) 

A A 

T = U 

The function, f($w), is the inverse of 6w(t^). In the program at this writing, the 
weight flow rate is a piecewise -linear function of t^, which permits a closed-form 
expression for f(6w). 
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3. 1.2 Constraints 

The arrival state, X^, is an implicit function of U, obtained through 
trajectory propagation from midcourse to target closest approach, which occurs 
at time t^. The constraint errors are defined explicitly as functions of desired 
(denoted by subscript d) and actual target (i. e. lunar) arrival conditions. 


\ = B * T d - B*T<X f ) 


0, = B'R _ - B*R(X,) 
2 d v r 


*3 * l fd - l f 


*4 = V hed - V he (X f> 


= V - V (X„,U) 
r 5 prd pr v f 


«/> 6 = m 


t = r , - r <XJ 
7 pd p f 


*8 ■ ‘d - ^ 


Equatorial miss -vector component 
Perpendicular miss -vector component 
Time of flight 
Hyperbolic excess speed 
Post-retro speed, periapsis maneuver 
Total correction fuel 
Radius of closest approach 
Inclination to target's equatorial plane 


Although MAESTRO allows targeting to input values of B- T^ and B- R^, it is more 

suitable (relative to orbiting mission objectives) to target to r ^ and i^. Convergence 

is surer when targeting to B*T^ and B*R^, however, so these are defined in terms 

of r , and i,. The specific definitions are found in appendices B, G and Hi. Nulling 
pd d 

of ^ and 

B«T, B-R, v T and r as functions of X„ can be found in the description of subroutine 
he p f 

BVE. We will now proceed to define v^ r (X^, U), i(X_^) and f(U). 


is a sufficient condition for nulling and 0^, The explicit definition of 
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The post-retro speed is defined by 


V = V °3 + 2 ^ r pd ' AV r 

which implies that the retro-velocity impulse, AV^ , is fired when the radius is 

r and opposite to the velocity vector there. The retro-velocity impulse is a 
pd 

function of the midcourse velocity correction (i. e. of U) since the retro engine is 

a solid of fixed impulse and U determines how much fuel remains at retro time. 

When v , is Ja/r ", and 4 = = 0 , the post-retro orbit will be circular, 

pra * pel a> 7 

Inclination is defined as that of the post-retro orbit relative to the target's 
equatorial plane, and therefore depends on the retro strategy adopted (see appendices 
D and E). If the angular momentum of the post-retro orbit is H, inclination is 
defined by 


i = cos 

where h is the polar component of H. 

O 

The total correction fuel, f(U), is the sum of the midcourse fuel and the trim 
fuel required to circularize the final orbit at the desired altitude. The midcourse fuel 
is computed from the rocket equation, assuming U to be the midcourse correction 
velocity impulse. The trim fuel is computed from the trim velocity, again using the 
rocket equation to provide the transformation from velocity to fuel weight. The trim 
velocity requirement is influenced by the retro strategy adopted and by the trim 
strategy. The normal strategy is to burn the retro at approach periapsis and to 
trim with two Hohmann impulses plus a plane-change, if necessary, at the larger- 
radius node on the target's equator. 
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3. 1.3 Guidance Laws 


Five guidance laws are considered. Each is designed to satisfy a different 
set of mission objectives, as seen from the following table. 


Guidance Laws 


Law 

Abbreviation 

Constraint Errors Nulled 

Minimum Midcourse Fuel 

MFG 

ih v tb z . |uf 

Fixed Time of Arrival 

FTA 

i!>2> 

Fixed Target Energy 

FTE 

V *2' *4 

Variable Target Energy 

VTE 

*!■ >V % 

Minimum Total Fuel 

MTF 

W*# ' 


* Minimized, rather than nulled 

The particular value of the MFG and MTF laws is that they achieve the final 
desired orbit with the least fuel. The FT A law is best from tracking visibility 
considerations, assuming that the arrival time was designed for good visibility. 

The FTE law is needed only if the arrival phase has been planned for a particular 
arrival energy and is inflexible to variations. The VTE law is useful in eliminating 
major trim maneuvers at the expense of midcourse and, usually, total correction 
fuel. 
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3.1.4 Pre-targeting 

The pre-targeting method to be described is used as a device to provide 
"first-guess" midcourse corrections to initiate precise differential-correction 
targeting procedures. In this role, the method eliminates the need for non-linear 
targeting measures such as control limiting and gradient re -computation and 
greatly reduces the number of required trajectory computations. A gross description 
of the method would be "patched-conic targeting with constrained end conditions. " 

Transfer Phase 



The transfer phase is treated as a single conic section relative to the 
central body (Earth). This phase begins at the midcourse correction position and 
time ends at the sphere of influence of the target body. The first terminal of this 
phase is the point of the midcourse correction on the uncorrected transfer trajectory — 
at radius R^, velocity and time t^. The second terminal is on the sphere of 
influence of the target at a radius R (referred to the central body) and time t . 

C* Lt 

The second terminal's radius vector is computed from the target's state (relative 
to the central body) at t , R and V 

Z t t 

R 2 ' E t + P 

where P is a vector from the center of the target body to the point of entry of the 
target's sphere of influence. The time, t , of arrival at the sphere of influence 

u 
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can be computed from 


where is the time of arrival at closest approach {specified a priori) and where T 

is the time required to travel on the approach hyperbola to (die point of closest approach, 

The arrival phase parameters, P and T, will be described later. An iterative solution 

of Lambert’s problem (i_e. "Find the conic section passing from R^ to R^ in time 

t - t ") is employed. This solution provides V , the velocity on the transfer conic 
2 1 ^ 1 

at the first terminal, and V the transfer conic's velocity at the second terminal. 

The midcourse correction impulse, AV, is computed from 

AV = V* - V, 

1 1 


and the target- relative velocity, V, at is 


V = V 2 - V 


where V is the target's velocity at t relative to the central body, 
t 2 


The first transfer conic is computed with r initially set to 66 , 000 seconds 
and P oriented in the Moon's orbital plane 10° earthward from the Moon's velocity 
direction. The target-relative velocity derived from this first conic solution is 
used to initiate the arrival phase calculations. The arrival phase calculations predict 
a new value of P to which the transfer phase is targeted. An iteration initiated in 
this way can be made to converge to a steady-state value of AV. 

— Transfer phase — 

P V 

Arrival phase 


This answer is interpreted as the impulsive-correction, patched-conic, fixed-time- 
of- arrival midcourse maneuver. 39 



Arrival Phase 

The arrival phase computations use the target-relative approach velocity 
and the desired end conditions to develop the point of entry onto the sphere of 
influence and the time from "patch" to closest approach, assuming the approach 
trajectory to be a target-centered hyperbola. The desired arrival conditions 
are specified values of radius of closest approach and inclination. We develop 
characteristics which the approach hyperbola must possess in order to satisfy 
the desired arrival conditions. We assume that the target-relative approach 
velocity vector defines the direction of the arrival asymptote of the approach 
hyperbola. Furthermore, we assume that the point of entry into the sphere of 
influence can be changed without changing the direction of the arrival asymptote 
or the target-relative energy. It is this last assumption's inaccuracy that prompts 
the iteration. 

The target-relative energy (Jacobian energy) is computed from the target- 
relative approach velocity, V. The energy of the approach orbit is corrected for 
perturbations which occur during the transfer phase but are not accounted for by 
the conic transfer trajectory. 

c 3 = V>V - 2 h/sl 

In the correction term, is the Moon's gravitational constant and a' is the 
semi-major axis of the Moon's orbit about the Earth. Note that the term (2it/r .) 

SOI 

is absent from the calculation. When the fully-perturbed trajectory is flown to 

lunar closest approach, the energy computed there is always weli -approximated 

by C of the above equation. Together with desired values for closest approach 
^ * 
distance, r and inclination, i^, the energy and velocity direction, V, can be 

used to compute the "desired" approach hyperbola and its corresponding point of 

entry into the sphere of influence, P. 
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The various parameters of the hyperbola are computed as follows. 



a 


-^■3 

- J C 3 + 2fi/r pd 

' ^VpV 111 

= -M/Cj 


b = (r pd - a) sin a 
e - J 1 + (b/a) 2 

. 2 

p = - a (b/a) 

n = /ju/|a|/|a| 

f = cos 1 ((p/66000“l)/e] 

M = TRMN(l,f, e) 

r = M/n 


hyperbolic excess speed 
closest approach speed 
half-angle between asymptotes 
semi -major axis 
asymptotic miss distance 
eccentricity 
semi-latus rectum 


mean motion 


true anomaly at patch distance 
mean anomaly at patch distance 
time from patch to closest approach 


We next calculate the point of entry into the sphere of influence in such a way 
that the hyperbola has the desired inclination, i We first define the arrival asymptote, 

^ A 

S, as the direction of the target-relative velocity at patch, V. Then, if K is a unit 
vector normal to the target's equatorial plane at time of closest approach, we can 
define vectors T and R normal to S as follows. 


T=SxK/SxK| 
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The miss -vector, B, lies in the orbital plane and in the plane of T and R at an 
angle, 0, measured from T toward R. 

B = T cos 0 + R sin 0 

The plane of the hyperbola is defined by a unit vector, H, in the direction of the 
angular momentum. 

H = B x. S = R cos 0 - T sin 0 

The condition that the orbit's inclination is i, is 

d 

H-K = cos i, = R*Kcos0, 
d 

which can be solved for 0 if j R - K j s jcos i^j. If |r-k| > |cos i^J , it means that 

i, cannot be attained for the S under consideration. In that case, the best that can 
d 

be done is 

cos 0 = 1 -sign(cos i^/R*K). 

The sign on sin 0 ( = ± J 1 - cos 2 0 ) can be chosen to make the miss-vector lie 
above or below the equatorial intersection in the miss-plane. Having now calculated 
B, we can form the vector, P, from the target's center toward the point of entry 
into the sphere of influence. 

P = Scos(f + o;) + Bsin{f+0') 

The computations above provide the point of entry of the sphere of influence and the 
time of passage to closest approach, which then can be used to establish a new aim 
point for the transfer trajectory. 
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Guidance Laws 


The solution provided by the above process Is the fixed time of arrival (FT A) 
guidance solution. The FT A guidance law constrains flight time, t f , to closest 
approach to be a specific value while satisfying the desired end conditions of radius 
of closest approach and inclination. Other guidance laws of interest relative to the 
RAE-B mission are: 

1 . minimum fuel 

2. fixed target energy 

3. variable target energy 

Each of the other guidance laws constrains radius of closest approach and inclination 
just as the FTA law does, but doesn't specifically constrain flight time. The minimum 
fuel law embraces the critical plane solution, the fixed target energy law constrains 
hyperbolic excess speed at the target and the variable target energy law constrains 
post-retro speed subject to a prescribed de-boost strategy. The solution 
for each of these laws, however, corresponds to a particular flight time, so flight 
time is used as the independent variable in seeking a solution for each law. 

The MFG law is targeted by means of a Newton-Raphson-type iteration with 
flight time as the independent variable. The iteration seeks to null the dot product 
of the difference of two successive midcourse correction impulses with the impulse 
itself. That is, it seeks to find the flight time for which the magnitude of the 
correction velocity doesn't change (i. e. is minimum). The condition for minimizing 
the magnitude of AV is 

= o. . 

dt. 


We define t to be the flight time for the n-th trial and AV to be the impulse for 
n n 

that trial. Then, approximately, 



AV - Av , 
n n-1 


t - t 
n n-1 
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and, more approximately, 


AV 


n+1 



<t , - t ) 
' n+1 n 



n 


We look for t such that the following equation holds. 



By making liberal use of the stated approximations, the preceding equation may 
be solved to render the time-step for the next iteration step. 


n+1 


- t - 

n 


- (t - t J 
n n-1 


AV 


n 


{AV -AV J 
n n-1 
2 — 


AV - AV 

“ n n-ll 


The process converges well for all cases tested, provided that the iteration for 
AV^(t^) converges well. The resultant maneuver changes very little with precise 
targeting, although the corresponding flight time may shift by an hour or two. 


The fixed time of arrival guidance law is pre-targeted to 1.015 times the 
desired flight time. This empirical factor tends to compensate for the difference 
between patched-conic and integrated-precise lunar transfer trajectory flight times. 


The pre-targeting process for the FTE and VTE guidance laws includes a 
standard Newton-Raphson iteration to null the third constraint by varying flight time. 
For the FTE law the third constraint error function is 





while for the VTE law the function is 


*3 = • / V+ pd - J c 3 + /* pd * 8v r 


M 



where 6 v is the velocity impulse imparted by burning the retro motor. There are 
usually two flight times for which = 0 for the VTE law. The iteration is constrained 

O 

to find the solution with positive slope — which is the solution corresponding to the 
longer flight time. If the minimum fuel solution is characterized by ^ >0 (the overbum 
situation), it will not be possible to find any VTE solutions. If this is the case, 
the precise targeting mode will probably not find any solutions either. 
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3.1.5 Targeting Methods 


In the neighborhood of U» variations in U are related to variations in jfiby 
5i/t = (7 ) fy U = G6 U 

The generalized Newton-Raphson method is used in the search for U* such that 
ib (U*) =0. If, at the k-th step of the iteration, we compute 


u k - Vi +5U 


" U k-1 - G k-1 


then will tend toward U*. The gradient, G = ib , is re-computed by finite 
differences (the Secant Method) at each of the first n iterations (n pre-set), then 
held fixed. To test convergence, the individual elements of ^ are tested against 
pre-set tolerances. 


The straightforward targeting method just described is used with the FTA, FTE 
and VTE laws, although the convergence of the VTE process is aided by certain 
modifications to G, 1 and ’boot-strapping" because 0 is an explicit function of 
U for this law. The MFG and MTF laws are treated as special cases, but are 
both formulated to use the Newton-Raphson method to obtain minimum-fuel 
solutions. In the case of the MFG law. 


U 


mm 




where C is the third column of G 


row to 


Ma 


and 


d iZ) «■ 


— 1 

k-1 * 


G, 1 is made invertible by adjoining a 


which is normal to them both; namely. 


dib 


1 


^2 


su au ' *" su " 'ju 

The MTF solution is also a constrained minimum. It is obtained by using flight 

time, T, as the independent variable in a scalar search for minimum A,. That is, 

6 

the FTA law is used to find a solution for each T, which renders A (T), i b^~\b n ~ 0. 

T is stepped away from the MFG solution time until ^(T ^) ^^(T ) £0 Q (T +1 ). A 

parabolic fit is then used to obtain T* for minimum i b„. The MTF solution is then 

o 


found by targeting FTA to T* . 
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3.1.6 Retro Maneuvers 


The retro or de-boost maneuver of the RAE-B mission is required to transfer 
from a lunar approach trajectory (which is approximately hyperbolic) to a near- 
circular capture orbit. The retro motor is a solid rocket whose firing direction 
and ignition time can be controlled, but whose impulse is fixed. The velocity 
imparted by firing the retro motor depends on the spacecraft's weight and, 
therefore^ on the midcourse fuel expended. This velocity can be more than 
sufficient to circularize the capture orbit at the desired radius (overburn), although 
for the RAE-B weights and specific impulses it is expected to be insufficient 
(underburn). The retro firing strategy (i.e. , specification of firing direction and 
ignition time) has an influence on the midcourse maneuver strategy as well as 
on trim strategy. We will proceed to discuss various retro strategies and their 
impact on midcourse and trim requirements. 

The objective of mission control for the RAE-B mission is stated here as achievement 

of a circular orbit of specified radius, r^, and inclination, i^. This objective, 

particularly the inclination part, can best be achieved with a midcourse maneuver. 

The midcourse maneuver should be targeted to a radius of closest approach of r^ 

and inclination of ij for underburns. For the overburn situation, several different 
d 

retro strategies can be defined to minimize total correction fuel requirements. Each 
retro strategy influences the midcourse targeting strategy. We now consider two 
midcourse targeting procedures which are applicable to overburns and involve the 
retro strategy. 

Variable Target Inclination 

The VTI procedure is derived in Appendix B. It amounts to defining desired miss 

vector components, B- T and B-R, in such a way that an out-of-plane retro-burn at 

peri apsis of the approach hyperbola can render a circular orbit of the desired radius 

and inclination. This procedure is applicable only for overburns. The sketch shows 

the retro velocity impulse, £V , in the plane of the approach periapsis velocity, V , 

r p 
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and the approach angular momentum vector, H. 


H 



Geometry of the VTT Procedure 

a 

The direction of AV^ is given by 

A A A 

AV t - Vp cos a + H sin a 


and the new angular momentum vector's direction is 

A ^ A A 

H = Hcos B - V sin B , 

P 

The angles, or and 0, are defined in the derivation of the VTI procedure as functions 

of r , d , and the desired final orbit radius. The desired miss parameters are 
3 v 

computed in such a w r ay that the approach orbit's track differs by £ from the desired 
orbit's track at r^. If the closest approach and inclination constraint errors are 
exactly zero, and if the retro-burn is perfect, the post- retro orbit will require no 
trim. 


Variable Approach Distance 

The VAD procedure, like the VTI, is applicable only to overburns. The basic idea 

is to define the desired miss vector magnitude in such a way that if the retro is 

fired appropriately in-plane at the desired circular orbit radius, r , the post-retro 

c 

orbit will be circular. It is assumed, first of all, that closest approach radius, r^, 

of the arrival hyperbola can be varied without significantly changing the arrival 

energy, C„, or the retro velocity impulse, d v. Assuming further that r is less 

3 p 

than the desired circular orbit radius, r , we can write 

c 


V (r ) = v R + v 9 
' c 7 r c Q 
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where v is the radial component of velocity and v^ is the tangential component 
The circular post retro velocity is: 


V (r ) = v 9 
c v c c 


where v 

c 



The required retro impulse, AV, is 


4V = V c - V = - v r R + (v c - v 0 ) 9 


2 2 ‘22 2 
and = v +(v - v J = v - 2 v v^+v 

r c u c 9 c 


C + 

3 r 


- 2 


{ T 


r + r 
c c 


= C + -2 -2- J C + — 

3r ^ r r » 3 r 

c c c p 


The only unknown in the above equation is r . We can solve for r as follows. 
J P P 

2 

2 U 2 _ 3 u i c i .2 

C 0 r +2, 

3 p 


" )[ 6v2 - C 3- V 4 -] 2-^| =» S 




_ J!_ 

C 3 


The periapsis constraint error is re-defined for this case as: 

PSI (7) = RP - PR. 

The desired miss vector magnitude is easily formulated from r . 

( ' i v v r \ P . 

tan” 1 P • °° - - P - ) 

v ** / 




The true anomaly at which the retro motor is fired is found from 



where the sign of 9 is determined by input: if IVTI is positive, 9 is negative 

A A 

and vice versa. Defining unit vectors I? and Q along periapsis position and 
velocity vectors, respectively, 

a A A 

R = P cos 9 + Q sin 0 

a a A 

9 = -P sin 9 + Q cos 9 

We define A V in terms of an out-of-plane angle, a, and a flight path angle, y, 
r 

If the periapsis constraint error is zero, a - it , and the maneuver will be in-plane. 

If 6v is too large to circularize at r for the actual r and C„ , the maneuver will 
& c p 3 

be out-of-plane. The direction out-of-plane is chosen to minimize the resultant 
inclination error. If 5v is too small, y will be chosen so that the radial component 
of post-retro velocity will be nulled, with the tangential component falling wherever 
it may as a result. 


smy = - 


cos a - 


( v c - v 0 


6v cosy 
r 


or a = IT if f cos Of [ >1 


AV f = Sv^ j(0 cos a + H sin aft cos y + R sinyj 


It is usually rather inexpensive to change closest approach distance by a few hundred 

kilometers. In such a change, target-relative energy is also rather invariant. If 

the final circular orbit radius is not very critical, it may be possible to change the 

targeted value of closest approach radius so that an in-plane periapsis retro maneuver 

opposite to the velocity vector can circularize the orbit. This discussion derives the 

value of r to aim for. 

P 


Circular Orbit at Arbitrary Radius 


Given and 5v, find peri-radius at which periapsis in-plane retro maneuvers 
render circular orbits. 


v c~ l v p " 6 v ) = ° = ~ V C s + 


+ 6 v =0 
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Finally, after some intermediate algebra, 

^ = [dv + V 2 5 v 2 - c 3 ] 2 

Approach Retro Strategy 

Following the final midcourse maneuver, residual approach errors are expected, 
due to pre-midcourse navigation errors and midcourse execution errors. Provided 
that the approach orbit is not expected to impact the Moon, quite a few retro firing 
options are available. Numerical methods are available within the approach analysis 
capability of MAESTRO for determining the retro controls to minimize trim 
requirements, but it is of value to examine and understand certain of the possibilities. 

It is believed nearly optimal from trim considerations to correct inclination errors 
as much as possible with the retro maneuver. The determination of optimal out- 
of-plane maneuvers is a complicated subject and is not covered here. A retro 
maneuver to correct residual inclination errors can be easily developed if the 
maneuver is executed at the node of the approach orbit on the lunar equator. An 
inclination-correcting maneuver at an arbitrary true anomaly can also be developed 
quite easily if certain assumptions are made on flight path angle of the maneuver, 
but will also not be elaborated here. 

Retro-maneuvers which are constrained to lie in the plane of the approach trajectory 
can be much more easily analyzed than are inclination-changing maneuvers. The 
following discussion treats in-plane circularization at an arbitrary radius using the 
retro only. 
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In-plane circularization at an arbitrary radius using the retro only 


Given: approach hyperbola (C , r or equivalently, C , h) 

op O 

retro impulse magnitude, <5v 
The velocity at radius r is 

A A A 


V = v r R + v Q e 


(9 is unit tangential vector) 


The circular velocity at r is 


V = v 9 
c c 


The velocity impulse required to circularize at r is 


AV = V - V = - v R + <v - v ) 9 
c c r ' c 9 



Let us examine the above equation in detail 


+ 3h «PV 

dr r^ ' 





^( r p 2 C 3 + - 2 JU r p ) (*^-) ~H 


At r = r^ , the slope is positive, because 


C +2u - ^ >0. 


P 3 


The slope goes to zero {a maximum of fjv ) at 

c 


r 


max 




+ 2 r 



A 

and remains negative (though ever less so) as r -4 . The limit of 6v c as 

r-toisC . If r C„< , the limiting fiv will be less than fiv (r ). Otherwise, 

3 p 3 4 c „ cp 

2 

6v (r ) is the lowest value of fiv . The curve of fiv vs. r has the general appearance 
c' p' c c 

shown in the following sketch. 



What does all this mean about retro-firing strategy? First of all, it means that if 

r Cg > , as is likely for RAE-B following the last midcourse correction, the 

’’cheapest" radius at which to circularize is periapsis. Then, if the available impulse, 

fiv, is greater than fiv (r ), it is possible to circularize the orbit with the retro alone 

c p 

at some radius, r, (r > r ). This situation may be called "an overburn at r ." If 

P - • P 

*v is also greater . than , an additional circular orbit radius may be found (r>r ). 

o mflA 

Of course, either of these circularization situations requires fiv^rfiv^r^^). If 

fiv is less than fiv (r ), the orbit cannot be circularized with the retro alone at any r, 

P 

but the "most circular" orbit is obtained by firing the retro at periapsis (unproven, 
but inferred). Even if the orbit can be circularized with retro-fire alone, the radius 
of the orbit may not be acceptable. Let us denote a desired circular orbit radius by 
r Let us also consider several cases which might arise. (C r > 

Case 1: r >r, 

P d . 

No circularization is possible at r^. 

- • 

la: fiv <fiv (r ) 

lLe 

No circularization is possible at r , but the maneuver at r will minimize the 

P P 

required trim. 

lb: fiv (r ) £ fiv 5 fiv (r ) (overburn at r ) 

C p C IUcIa p 


A circular orbit can be obtained at r , where fiv (r ) = fiv. Now, r can be 

c c c c 
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compared with - it may be an acceptable radius. At r > r c> 6v represents 

an underburn at r, while at r in the range r s r < r , 6 v cau provide an overborn 

P c 

at r. In the latter case, r is the apoapsis radius of the post-retro orbit and it 
may be possible to lower post- retro periapsis to r^; a minimum-trim strategy. 

It is not possible to lower post- retro periapsis below that which would result 
from overburning at r . These recent statements assume that the post-retro 
velocity lies entirely in the tangential direction, i.e. , along 0 . It is also possible 
to give the post- retro velocity an out-of-plane angle and/or a flight path angle, 
thereby using fiv less efficiently and changing the post-retro orbit's shape and 
orientation. 

lc: 5v > fiv (r ) 
c max 7 

► 

Not likely to occur. 


Case 2: r <. r. 


Three regions can be identified. The first is when §v is greater than necessary to 

circularize at r . In this case the retro can circularize the orbit at r (r <r < r ). 

p c p c d 7 

If 6v</5v (r ), the "most circular" place to fire the retro is at r . The best place 
c v p 7 p 

to fire from the criterion of minimizing trim to r^ depends on the opposite apsidal 

radius, r , which is achieved by firing the retro at r . 

a p 

2at r >r, 
a d 

The best place to fire to minimize trim is r since the change of r^ with r has 

a slope of about 7 for likely RAE-B ranges of C , r and <5v. This slope implies 

3 p 

a trim fuel penalty for firing before or after approach periapsis. 

2b: r <r, 
a d 


Fire the retro at r such that r = r,. This will minimize the trim. 

a d 

See section 3.3 for more information on retro strategy. 
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3.1.7 Lunar Orbit Trim 


Once the retro motor has fired, the spacecraft is in lunar orbit. Assuming this 
orbit to be elliptical with energy C^ e and angular momentum h^, we will proceed 
to define the trim sequence and the trim fuel cost. The in-plane trim is accomplished 
with a 2-impulse Hohmann transfer to circular at the desired radius, r^. 



The 


periapsis radius, r , and apoapsis radius, r , can be defined from C 

p a 


3e 



semi-latus rectum 


and 


e 

r 

P 

r 

a 

a 


= t|i +-Si<£ 

M 

^ 

1 + e 

. £ 

1 - e 

, r P +r »- s 

2 



'^eccentricity 
periapsis radius 
apoapsis radius 
semi-major axis 



Case 1: r > r, 
a d 


The semi-major axis, a., and energy, C , of the intermediate orbit must be 

1 ol 

r a + r d 
a = — 
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and 


JL- = _ -llL 


3i 


a, 

1 


r + r . 
a d 


The velocity at r before the first impulse is 
a 


V = lie + 2lL. = J ~ r P 

a v 3e r 1 v /r . , 


r (r +r ) 
a' a p' 


and after the impulse, 


a 


C 3i + f- 
a 




iM 


r (r + r.) 
a' a d' 


so that the first maneuver requires an impulsive velocity 


v = ^ I - J 

1 y r r +r . y 

a ad 


r + r 
a p 


The velocity on the intermediate orbit at r , is 

d 


C 3i + r 


l tL 


= 

’ r .(r + 


,(r +r J 
d' a d' 


while the desired circular velocity is 

v* = 
d “ r . 

d 


so that the second impulse is 



and the total in-plane trim velocity is 
6v = 6v 1 +6v g . 


Case 2: r^cr^. In this case, it is slightly cheaper 
By a derivation similar to that for case 1, we obtain 


to transfer first from r to r 

P d 






It should be noted that the two cases can be computed with the same equations by 
interchanging roles between r and r . 

cL P 


Inclination Change 

The inclination is changed by a third impulse executed at the node of the orbit on 
the equator. The impulse is applied on the intermediate orbit of the in-plane 
adjustment at the nodal crossing of larger radius. Although this strategy is not 
usually optimal, it is convenient, it leaves the node invariant, it always permits 
a solution and it separates out the in-plane and out-of-plane trim costs. Let the 
larger nodal radius vector (of either the pre-trim or intermediate orbit, . actually) 
be R and its corresponding velocity be V . If the inclination is to be changed by 
fii, the rotated post- impulse velocity, will be 


V 


+ 


and 6 V 

O 


cosfiiV + (1-cos 6 i) R*V R + sin 5 i R x V 

+ -* 

= V - V 


= (cos gi - 1) V + (1-cos fii) R • V R + sin 6i R x V 

2 r 2 2 1 r - - A - 2l 

J [v * v - (R-v ) J 


fiv, 


- 1 (1 - 


2 2 
cosfii) + sin 6iJ 


. A - 1 2 

= 2 (1-cos 6i) |R x V j 


= 4 sin' 


. 2 / ai \ 
m 


6v = 2 sin 

O 


M. 


R x V 
|r x V 


To roughly estimate the cost of correcting inclination, assume 

A*,. _ 

and sin ~ — = — - — . 


| H x v"| = 


6v 


« = 22.94 (JS/deg) fii 0 


( 
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The fuel cost is computed with the rocket equation 


6v = 6v x +6 v 2 +6v 3 



The weight before trim is w^ and the characteristic velocity 
motor is c. 


(gl ) of the trim 
sp' 
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3.2 Operational Strategy 

The capability of MAESTRO'S midcourse guidance package enables the Flight 
Director to design maneuvers enhancing the probability of mission success. 

This section is intended to describe how one might best use the program to 
answer the midcourse correction questions. 

3.2.1 Mission Success Criteria 

The primary criterion for mission success (from a guidance point of view) is 
whether enough correction fuel is available to place the spacecraft in a circular 
lunar orbit at the right altitude and inclination. The solution is constrained to a 
specific range of allowable angles between the Sun-direction and the spin-axis 
direction at any maneuver time and is also constrained to insure tracking vis- 
ibility everywhere in the vicinity of the maneuvers. 



3.2.2 Operational Procedures 

The procedures for real-time selection and determination of RAE-B's midcourse 
correction maneuvers will now be described. Figure 3 . 1 illustrates the order in 
which the program’s various operational modes can be utilized to select a suitable 
maneuver. The first step shown in the figure is propagation of the anchor vector 
to closest approach to the Moon and determination of trim requirements. This can 
be done by invoking the trajectory propagation (MODE = 0) and approach analysis 
(MODE = 2) modes. Alternatively, this step can be taken by running the midcourse 
analysis (MODE = 3) mode, since the first printout from the midcourse analysis . 
describes arrival conditions and approximate trim requirements for the no-mid- 
course trajectory. The second step is to run a minimum-fuel guidance (MFG) scan 
for a range of allowable midcourse ignition times. This scan will determine total 
correction fuel requirements, spin-axis/Sun angle and visibility for midcourse and 
retro maneuvers for the MFG law. The maneuver will probably be selected from 
the FT A -law scan, presuming that arrival time was selected for best visibility, 
but the MFG scan is useful in ascertaining primary mission success attainability. 

The primary criterion for mission success is required- vs. -available correction fuel, 
and if the mission objectives cannot be achieved under the MFG law, there is little 
chance they could be achieved under a different law. The midcourse analysis 
capability does not include a facility for studying alternate mission possibilities, 
although marginal opportunities may be enhanced by varying targeting and retro- 
strategy parameters which are specified by input. Assuming that the available 
correction fuel supply is more than sufficient to achieve the desired orbit under 
the MFG law in an attainable range of midcourse ignition times, the next step 
(3a in the figure) can be taken. 

Scanning runs under the fixed-time- of- arrival (FTA) and variable -target-energy 
(VTE) laws will obtain the penalties in fuel and penalties or advantages in Sun angle 
and mane uver visibility of imposing arrival constraints. If the fuel cost is acceptable 
and the Sun angle constraint is not violated, the FTA law would seem to render the 
most likely candidate maneuver from visibility considerations, assuming the correction 

is targeted to the nominal arrival time. If a fairly-large midcohrse maneuver is 

00 



Operational Strategy Diagram 


Figure 3.1 












required for the MFG and FTA solutions, the VTE solution can be obtained without 
much penalty in fuel or visibility* This solution has the advantage of circularizing 
the lunar orbit with the retro-burn alone if the midcourse and retro burns are accurate 
enough. By eliminating a major trim maneuver, greater reliability should be achieved. 
There seems to be no advantage to scanning the requirements for the fixed-target- 
energy (FTE) law and very little advantage in the minimum-total-fuel (MTF) law 
for the real-time midcourse analysis of the RAE-B mission. 

The choice of a suitable midcourse correction maneuver must be made by the Flight 
Director, who makes his decision based on the information provided him. If only 
a small correction is required, he may somewhat arbitrarily choose a late maneuver 
and perform more analysis in the ample time left before the maneuver is to be 
executed. If a marginally -large correction is required, it will usually be advantageous 
to execute the first midcourse correction as early as possible. In this case, extensive 
analysis may have to coincide with or follow transmission of the maneuver command. 

The selected maneuver should be re-targeted with as much precision as is available 
in trajectory propagation to ascertain the specific values of burn time and attitude 
to be commanded. Once a suitable maneuver has been selected, the midcourse 
verification (MODE = 5) and approach analysis (MODE = 1) modes may be invoked to 
study the solution. Further midcourse guidance analysis (MODE = 3) capabilities are 
available, also, for illuminating characteristics of the solution. By scanning flight 
times (FTA) as well as midcourse ignition times, the Flight Director can gain an 
understanding of how the selected maneuver fits into the set of possible maneuvers 
and how variations in flight time or ignition time influence fuel requirements and Sun 
angle behavior. The "fixed- attitude" scan enables the effects of execution errors to 
be determined. Under this option, the program automatically varies midcourse 
correction velocity and/ or right ascension and declination of the thrusting axis over 
specified ranges at each midcourse ignition time. The output includes radius of 
closest approach, inclination, flight time and total correction fuel, enabling a 
determination of execution error effects on these quantities. This option also enables 
an understanding of penalties for incomplete attitude control which may drive the spin- 
axis to an attitude near-but-not-at the attitude required for a perfect correction. 
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The Monte Carlo (MODE = 4) capability can be utilized to assure the Flight Director 
that if he can believe the input tracking and execution error statistics, the success 
of the mission is assured. If the Monte Carlo run indicates low success probability, 
the assumed error statistics must be reduced. 



3.2.3 Run Setups 

A great many of the program option keys can be set up before the flight and run in 
blocks according to the operational procedure being followed. The launch epoch, anchor 
epoch and anchor vector must be provided in any case, as must also the attitude, 
weight and thrusting information. As for the remaining inputs, many of the inputs 
peculiar to midcourse guidance analysis are invariant with chosen procedure or 
guidance law and will be provided in Block Data. Library files containing the inputs 
peculiar to each procedure step should be set up before launch to enable rapid turn 
around of midcourse guidance information. It is the purpose of this paragraph to 
define such files. 

File 1: Epochs, anchor vector, attitude, weights, thrust tables 

File 2: Standard midcourse settings 

File 3: Settings for a step of the midcourse procedure 

Figure 3.2 Input data stack 

Standard settings for midcourse analysis 

Table 3.1. contains those settings for midcourse analysis which would be invariant 
between procedure steps. 


Table 3.1 

Standard Midcourse Settings 


410-419 

OBSLON 

*** 

Tracker longitudes (deg) 

420 

PSID{1) 

2838. 

Radius of closest approach for targeting (km) 

421 

PSID(2) 

121. 

Selenographic inclination for targeting (km) 

422 

PSID(3) 

414000. 

Nominal flight time from launch (sec) 

423 

PSID(4) 

.62 

Nominal at lunar arrival (Ion 2 / sec 2 ) 

424 

PSID(5) 

0 . 

Desired circular excess speed (km/sec) 

435 

SIGATM 

.5 

Expected midcourse attitude error (deg) 

436 

SIGDVM 

.01 

Expected midcourse proportional error (fraction) 
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440 

TMC2 

259200. 

Second midcourse time from launch (sec) 

444 

RO 

2838. 

Desired lunar orbit radius (km) 

446 

CIO 

121. 

Desired selenocentric orbit radius (deg) 

449 

TRINC 

0. 

Inclination trim tolerance (deg) 

475 

TRUE 

20. 

True anomaly range for retro firing (deg) 

480-89 

ORSLAT 

*** 

Tracker latitudes (deg) 

490 

TOL(l) 

10. 

Tolerance on B* T (km) 

491 

TOL<2) 

10. 

Tolerance on B*R (km) 

492 

TOL(3) 

10. 

Tolerance on flight time (sec) 

493 

TOL(4) 

.0001 

Tolerance on hyperbolic excess speed (km/sec) 

494 

TOL(5) 

.0001 

Tolerance on circular excess speed (km/sec) 

495 

TOL(6) * 

.02 

Tolerance on total fuel (kg) 

496 

TOL(7) 

5. 

Tolerance on closest approach radius (km) 

497 

TOL(8) 

.2 

Tolerance on approach inclination (deg) 

1036 

KMETH 

6 

Trajectory computation method for analysis 

1044 

MODE 

3 

Midcourse analysis mode key 

1050 

MCOUT 

0 

Extra output key for midcourse analysis 

1058 

KOUT9 

9 

Auxiliary unit number (logical for scope) 

1061 

MCUNIT 

11 

Auxiliary unit number for temporary storage 

1062 

IBTR 

2 

Selector for targeting to radius and inclination 

1064 

NGROPT 

1 

Number of trials to re-compute the gradient 

1065 

NT 

10 

Maximum trials for targeting . 

1066 

JET 

1 

Preliminary targeting key 

1067 

MCLIM 

100 

Limit on step-size for targeting 

1070 

IPROB 

0 

Probability for linear propagation (%) 

1071 

IBURN 

6 

Midcourse burn method key 

1075 

KMETHP 

6 

Initial integration method for targeting 

1077 

KTF 

0 

Scan procedure key 

1078 

IVTI 

2 

Variable target approach key (set for in-plane) 

1080 

NORMIN 

0 

Retro maneuver optimization key 

1093 

KPLOT 

0 

Plot key for scanning modes 
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Minimum Fuel Guidance Scan 


The first run should be a broad scan of MFG requirements near anchor epoch to 
establish the general scope of midcourse and retro requirements. Table 3.2 
presents the parameters which must be considered for this step, together with 
appropriate values for them. 


Table 3.2 

Minimum Fuel Scan Step Inputs 


434 

DELTMC 

18000. 

Step-size for scanning ignition times (sec) 

478 

TMC 

18000. 

First midcourse ignition time (sec) 

1051 

JMC 

9 

Number of ignition times to consider 

1063 

KGLAW 

1 

Guidance law indicator (MFG) 


The results of this scanning run will be total correction fuel, spin-axis/Sun angle and 
tracking station visibility versus midcourse ignition time for the minimum mid course 
fuel solution. Answers will be computed with the multi-conic trajectory propagator 
and other simplifying assumptions, but will indicate the midcourse guidance requirements 
and choices. The first trial midcourse time of 5 hours may be too early to be achievable, 
considering anchor epoch and attitude -change and attitude determination time requirements. 
This run will hopefully reveal one or more ignition time regions in which a properly 
executed midcourse correction bum will render a good trajectory for total achievement 
of mission objectives. Violation of Sun angle constraints can usually be averted and 
poor visibility can be corrected or even accepted if necessary, but total fuel requirements 
are really the indicator of mission success achievability. If every ignition time of the 
scan required (for the MFG law) more than the available fuel, an ^valuation of alternate 
mission possibilities should be immediately undertaken. 

The standard setting for overburn retro-strategy (IVTT, location 1078) is 2. This setting 
obtains the in-plane pre-pericynthion firing strategy which circularizes the orbit with 
the retro-burn at the final desired radius. Other values for retro Sun angle and 
visibility can be obtained with settings of -2. , 1. and -1. Required midcourse fuel 
will vary only slightly with changes in the setting of IVTI and there is no difference in 

the solution unless the overburn situation is encountered (midcourse correction fuel 
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requirement is greater than about 16 kilograms). It is not known at this time which 
of the overburn retro firing strategies is preferable. 

Other Guidance Law Scans 

If the MFG scan indicates satisfactory fuel requirements for a range of ignition times, 
that range should then be scanned for FTA and VTE solutions. The fixed time of arrival 
run should be targeted to the nominal flight time to satisfy visibility constraints. The 
VTE law should be targeted only in the range for which the MFG solution is an underburn 
(6v less than about 16 kg). The following tables are set up for a case in which the 
MFG solution has a small fiv over the entire range of ignition times. The simplifying 
assumptions of the MFG scan are retained through this step. The FTA and VTE cases can 
be stacked. 


Table 3.3a 

Fixed Time of Arrival Scan Step Inputs 


422 

PSID(3) 

414000. 

Desired flight time (seconds from launch) 

434 

DELTMC 

18000. 

Step-size for scanning ignition times (sec) 

478 

TMC 

18000. 

First midcourse ignition time (sec) 

1051 

JMC 

9 

Number of ignition times to consider 

1063 

KGLAW 

2 

Guidance law indicator (FTA) 




Table 

.3.3b 


Variable 

Target Energy Scan Step Inputs 

424 

PSID(5) 

0 . 

Desired circular excess speed (km/sec) 

434 

DELTMC 

18000. 

Step-size for scanning ignition time, (sec) 

478 

TMC 

18000. 

First midcourse ignition time (sec) 

1051 

JMC 

9 

Number of ignition times to consider 

1063 

KGLAW 

4 

Guidance law indicator (VTE) 
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Verification of Candidate Maneuver 


After the Flight Director has chosen a candidate ignition time range, guidance law 
and retro firing strategy from the MFG, FTA and VTE scanning steps, he should 
scan ignition times in a smaller localized range and then re-target the particular 
solution without the simplifying assumptions used in the broader scan. 

This run defines the ignition time, burn duration and spin-axis attitude for commanding 
the maneuver. 


Table 3.4 

Verification Targeting Step Inputs 


422 

PSID(3) ’ 

414000.* 

Desired flight time (sec) 

434 

DELTMC 

3600. 

Step-size for scanning ignition times (sec) 

478 

TMC 

50400.* 

First midcourse ignition time (sec) 

492 

TOL(3) 

1 . 

Flight time tolerance (sec) 

496 

TOL(7) 

1 . 

Closest approach tolerance (km) 

497 

TOL(8) 

.01 

Inclination tolerance (deg) 

1014 

KINT 

5 

Numerical integration scheme key 

1029 

KOBL 

1 

Lunar oblateness key 

1036 

KMETH 

4 

Trajectory propagation method 

1051 

JMC 

3 

Number of ignition times to consider 

1063 

KGLAW 

2* 

Guidance law indicator 

1070 

KPROB 

99 

Probability level for linear error propagation (%) 

1071 

IBURN 

1 

Midcourse burn method 


* assumed characteristics of the selected maneuver 


The linear error propagation should really be run with the multiconic trajectory 
generator rather than in this step, but should not be run in a broad scanning step 
either. The selected maneuver should be used to initiate a run in the midcourse 
verification and approach analysis modes. 
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Flight Time Scan 


The flight time scan step will illuminate the nature of the solution relative to 
neighboring solutions. It would be better to run this scan before selecting the 
candidate solution if there is time to study the results before having to supply commands. 
An auxiliary data tape will be written as a step in generation of hard-copy plots of 
the scan results. 


Table 3.5 


Flight Time Scan Step Inputs 


422 

PSID(3) 

396000.* 

Initial desired flight time (sec) 

434 

DELTMC 

3600. 

Step-size for scanning ignition times (sec) 

478 

TMC 

50400.* 

First midcourse time (sec) 

1050 

MCOUT 

1 

Extra output key 

1051 

JMC 

3 

Number of ignition times to consider 

1063 

KGLAW 

2 

Guidance law indicator (FTA) 

1077 

KTF 

10 

Number of one-hour flight time steps 

1093 

NPLOT 

21 

Plot indicator and tape unit number 


* assumed characteristics of the selected maneuver 

h h 

This case will scan flight times from 110 to 120 in one-hour steps while varying 

h h 

ignition times between 14 and 16 in one -hour steps. Submittal of the plot-run 
should follow examination of the data. 
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Execution Variation Scan 


The execution variation scan can provide bounds on performance by perturbing 
the midcourse maneuver about its commanded value. The varied parameters 
(in order of most- rapid variation) are: 1. declination, 2. right ascensipn, 3. velocity 
impulse and 4. ignition time. 


Table 3.6 

Execution Variation Scan Step Inputs 


47 

RAI 

50.68* 

48 

DECI 

43.34* 

426 

PSID(7) 

.03807* 

434 

DELTMC 

36. 

474 

CONE 

1 . 

478 

TMC 

53964.* 

479 

DINK ' 

.001 

1041 

JRA 

3 

1042 

JDEC 

3 

1051 

JMC 

3 

1077 

KTF 

-3 

1080 

NORMIN 

2 


Desired right ascension of spin-axis (deg) 
Desired declination of spin-axis (deg) 

Initi al velocity impulse (km/sec) 

Ignition time scan-step (sec) 

Attitude scan-step (deg) 

Initial ignition time (sec) 

Velocity impulse step-size (km/sec) 
Number of right ascension steps 
Number of declination steps 
Number of ignition times to scan 
Number of velocity impulse steps (negative) 
Inclination-trim option key 


The variations shown in this table (.01 , 1°, 1 m /sec) are set to convenient values 
for interpolation, but could be set to 3o values to examine worst-case performance. 
After examination of the results of this run, the Monte Carlo option can be invoked 
to investigate statistically the effects of navigation errors as well as errors in the 
velocity impulse and attitude. 
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3.3 Experimental Results 

MAESTRO'S midcourse guidance capability has been tested relative to two launch 
opportunities (April 14, 1973 and June 10, 1973) using simulated anchor epochs and 
anchor vectors. The results of these tests are presented here to illustrate trends 
and are not to be regarded as quantitative results. . . 

3.3„1 Flight/Ignition Time Scans 

The option for scanning flight times with the fixed time of arrival (FT A) guidance law 
was used to generate midcourse maneuver requirements at each of several midcourse 
ignition times. Each maneuver is targeted to the same radius of closest approach- - 
and inclination. 

Figure 3.2 shows total fuel (i.e., midcourse plus trim) requirements, for the nominal 

April 14 trajectory. This figure shows 9 v-shaped curves. The uppermost of these 

curves describes total fuel for a midcourse maneuver executed two hours after epoch, 

with the other curves spaced at two-hour intervals thereafter. These curves focus 

h 

toward a point at 115.3 flight time and 6. 3 kg total fuel. This total fuel requirement 
is used entirely for trim, since no midcourse maneuver is needed for any candidate 
maneuver time when guiding FTA to the nominal arrival time. The minimum fuel 
guidance (MFG) law solutions for these same maneuver ignition times are all 
characterized by negligible midcourse maneuvers, 115. 3^ flight times and 6.3 kg total 
fuel. The figure shows the total fuel penalty for midcourse targeting to flight times 
. away from the nominal. 

Figure 3.3 shows the total fuel requirements for the April 14 opportunity when injection 
velocity is 10 m /s "hotter" than the nominal. The two-hour correction is represented 
by the curve with the lowest minimum. The minimum total fuel for this ignition time 
is about 8.4 kg and the minimum fuel increases with ignition time to about 13. 8 kg 
for the 18-hour correction. It may also be observed that the minima occur at flight 
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Figure 3 . 3 
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times slightly shorter than the nominal, a fact which correlates well with the fact 
of a "hot" trajectory. 

Figure 3.4 displays the same data as figure 3.3, although the two figures have vastly 
different appearances. Figure 3.4 is a contour plot which presents contours of equal 
total fuel in the plane of flight time and midcourse ignition time. The contour 
of least total fuel shown in the figure is the one at the bottom. This is the contour 
of 10 kg total fuel. Its irregularity is due to granularity of the data and to shortcomings 
of the contouring and plotting algorithms. The contour lines are separated by 1 kg 
in total fuel. The region of allowable flight times and midcourse ignition times would 
be bounded by the contour of maximum-allowable total fuel. 

Figure 3.5 shows still another way of presenting the data of figure 3.3 . In this figure,- 
total fuel is plotted against midcourse ignition time. Each curve represents the fuel 
requirements for a particular flight time, targeted with the FTA law. The uppermost 
curve is for 120 flight time, the next highest for 110 . The lowest two curves are 
for 114 a and 115 11 . These lowest curves indicate the total fuel penalty for delaying the 
minimum-fuel midcourse maneuver. 

Figure 3.6 displays the required spin-axis orientations for midcourse maneuvers on 
the "hot" April 14 trajectory. Each curve is for a single midcourse maneuver time and is 
generated by connecting points of different flight times. The obviously out-of-place 
lines from the uppermost curve are plotting errors and should be ignored. The upper- 
most curve is for the two-hour correction. The curves begin at the upper— right corner 
of the plot, where flight time is 110 h , and proceed toward the lower-left, where flight 
time is 120 . Curvature of the curves is exaggerated, due to plotting declination 
versus right ascension, but the path on the celestial sphere still differs somewhat 
from a great circle. The dot at about {20 , 10 ) is the direction of the Sun at midcourse 
time (it moves only about . 5° over the 16-hour span of midcourse times). The line 
drawn nearly perpendicular to the family of curves shows how midcourse spin-axis 
orientation varies with midcourse time for a 115 FTA correction (which is nearly 
the MFG correction for this trajectory) . 
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Figure 3.7 shows spin-axis/Sun angle, plotted against flight time, again for each 
of a span of midcourse ignition times. The figure shows the distance of the points of the 
curves of figure 3. 6 from the Sun-direction in that figure. The uppermost curve 
is for the 2-hour correction maneuver. 


Figure 3.8 shows contours of constant spin-axis/Sun angle in the plane of flight time 
and midcourse ignition time. This is another way of presenting the data shown 
in figure 3.7 . The lowest contour shown is the small v-shaped one at the top of the 
figure. It represents a spin-axis/Sun angle of 20°. Successive contours are separated 
by 5°. 


Figure 3.9 presents the spin-axis/Sun angle information of the previous two figures, 

but plots it against midcourse time in curves of constant flight time. The uppermost 

h h 

curve is for a 110 flight time and the next-highest for a 111 flight time. The 

curves get generally lower to the 115 flight time curve and then climb again. The 

lowness of the 115 curve shows that the minimum fuel correction has about the smallest 


spin ~ cuXi s /Sun angle and especially so fox* 1 a vO r midcourse ignition times (for this case). 


Figure 3. lOis the same type of plot as figure 3.6, but presents spin-axis orientation 

requirements for the April 14 "cold" trajectory whose injection velocity is 10 m /s below 

nominal. The curves are almost reflections of the curves of figure 3.6 about the 

0°-declination axis. The orientations of short flight times appear at the lower-right 

part of the figure and the curves focus toward the upper left as flight time increases. 

h 

The lowest curve is for the 2 correction time. 


Figure 3. IT is plot of total fuel versus midcourse fuel for the April 14 "hot" trajectory. 

Each sloping v-shaped curve is generated by scanning flight times for a particular 

ll 

midcourse execution time. The upper-left point of each cui’ve is for 110 flight time 
and the right-most point is for 120* 1 . The skinniest and longest curve in the center of 
the bunch is for the 2 midcourse time. The plot shows that the slope of total fuel 
to midcourse fuel changes on either side of the minimum fuel solution. It also shows 
that the minimum total fuel solution is very nearly the same as the minimum midcourse 
fuel solution. Other test results (see Progress Report for June, 1972), show that these 
curves focus into a straight line at the overburn point which begins at about 16 kg of 

midcourse fuel expenditure . If a circular i zing-retro overburn procedure such as 
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variable target inclination is used, the slope of this line is 45 , otherwise the slope 
o 

is greater than 45 . 

Figure 3.12 was generated for a "hot" June 10 trajectory. It presents trim fuel 
versus midcourse correction fuel for that trajectory. This figure presents the 
same type of information as figure 3.1 X. Each curve is generated by scanning flight 
times for a particular midcourse maneuver time. The curve containing the upper- 
leftmost point is for a 2 correction. Each curve focuses to zero-trim at about 16 kg 
of midcourse fuel expenditure and then bounces upward along a focused beam. If 
a circularizing-retro overburn procedure had been assumed, the post-retro trim 
fuel would have been zero and the curves would have followed the lower edge of the 
plot after about 16 kg of midcourse fuel expenditure. 

Figure 3.13 presents the same data as figure 3.12., but plots trim fuel against flight 
time. By extrapolating each of these curves to the flight time of zero trim require- 
ment, it is possible to obtain the flight time of the variable target energy (VTE) solution 
for each midcourse execution time. 

Figure 3.14 is a contour plot of trim fuel in the plane of flight time and midcourse time. 

The minimum contours shown are the third and fourth lines down from the upper right- 
hand corner of the plot. These lines of . 5 kg trim fuel bound a valley whose minimum 
is the zero-trim or VTE solution for this case. This contour plot also shows a maximum 
trim region in the lower left-hand corner, should a maximum be of any interest. 

Figures 3.15 and 3.16 present the relationships between midcourse fuel and velocity 
(delta-v) and between burn time and midcourse fuel expenditure. These curves are 
based on initial spacecraft weight of 333.39 kg, an I of 226, an initial weight flow 
rate of .02404 kg/sec and a weight flow rate at 900 seconds of .01202 kg/sec. One might 
use these curves as quick-and-dirty cross references. For instance, a 100 m/s correction 
could be translated into about 14.7 kg of expended fuel and then into a midcourse burn 
duration of about 13 minutes. 
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Figure 3.12 
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Figure 3.13 
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Figure 3.14 
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3.3.2 ” Fixed-Attitude 11 Scans 


The ’'fixed- attitude guidance” capability is simply an automated procedure for varying 
the midcourse command parameters: 

1. ignition time 

2. burn duration (A V) 

3. right ascension of the spin axis 

4. declination of the spin axis. 

Two situations have been investigated using this capability. The first of these is that 
the execution of the midcourse maneuver is imperfect relative to that determined by 
midcourse targeting. The automated procedure can be made to compute and display 
the consequences of such imperfections. The second situation which has been investi- 
gated is that the attitude control system is inoperative. The question of interest is then 
of what can be done with control of only ignition time and burn duration. 

Figures 3. 17 and 3. 18. show total fuel (midcourse plus trim) for variations of ±6 
minutes in ignition time, ±1 m/s in velocity impulse and ±1 degree in each angle relative 
to a 2-hour MFG-targeted correction. The fuel values are in kilograms and primarily 
reflect trim variations, since the midcourse fuel is nearly constant at about 6 kilograms. 
In computing the data shown in Figure 3. 1.7., the retro motor was fired in-plane at 
periapsis and inclination errors were not trimmed cut. For Figure 3, 18, the firing 
time and retro attitude were optimized to minimize trim velocity while trimming out 
inclination errors. Inclination of the post-retro, pre-trim orbit was printed out on 
the run and varied only occasionally and slightly from the desired orbital inclination of 
116. 5 . This result confirms the suspicion that it is nearly optimal to remove inclination 
errors with the retro burn rather than even partly with trim maneuvers. Figure 3.19 
shows the variations in radius of closest approach and inclination which result from 
command errors of the 2-hour MFG correction. An initially-surprising result shown 
by this figure is that attitude errors in right ascension and declination map into the same 
straight line on the figure. Appendix D explains why this result is generally true for 
attitude variations from a minimum fuel correction. Each arrow in the figure is the result 
of attitude variations at a particular ignition time and burn duration. The central arrow 
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Figure 3.17 
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Figure 3.18 
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is for the nominal ignition time and burn duration, while the extreme arrows are for 
±1 m/s deviations in velocity impulse from the nominal of 39.07 m/s. The other 
arrows are for ignition time variations of ±6 minutes from the nominal, holding the 
velocity impulse constant at its nominal value. An earlier-than-nominal ignition time 
is seen to have the same effect as a "hot' r burn and vice versa. 

Figures 3.2 0 3.21 , 3.22 and 3.2 3 describe the end conditions which could be 

achieved with the midcourse spin-axis attitude fixed at that of the injection velocity 
vector. The trajectory is the April 14 nominal 115-hour trajectory which requires no 
midcourse correction. The figures show separate curves for each midcourse ignition 
time under consideration. These are from 2 hours to 18 hours after launch. The top 
curve in each figure corresponds to the 2-hour correction. The abscissa in each 
figure is midcourse velocity magnitude in meters per second. Figure 3.2.0 shows radius 
of closest approach, figure 3.2 1 shows inclination and figure 3.2 2. shows flight time 
to closest approach. Figure 3,2.3 shows total correction fuel (midcotirse plus trim) 
for in-plane trim of the orbit to circular. The calculations used in generating this 
figure assume that attitude control becomes operable prior to the retro maneuver and 
that the retro maneuver occurs at periapsis and opposite to the velocity there. It may 
be noted that the variations are nearly all due to trim fuel requirements since the curve 
of midcourse fuel vs. midcourse velocity is very nearly a straight line through the origin. 

3.3.3 Retro and Trim Results 

Very few presentable results have been obtained concerning the retro maneuver. Solutions 
for retro direction and firing time are determined either by numerical optimization to 
minimize trim fuel or by an assigned retro strategy such as 

1. firing the retro at perilune to give an impulse opposite to the 
velocity there. 

2. firing the retro at perilune to impart an out-of-plane impulse required 
by the variable target inclination procedure, or 

3. firing the retro in-plane to circularize the orbit at the desired final 
radius in accordance with the variable approach guidance procedure. 

Figure 3. 24 is a contour plot which shows the path of the retro direction during the 
numerical optimization procedure. 
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Figure 3.2 3 
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Figure 3.2 



The numerical procedure for obtaining retro direction and firing time includes, 
for each candidate retro direction, a minimization of trim fuel by varying firing true 
anomaly. Figure 3.25 is a plot of trim velocity vs. firing true anomaly for two 
neighboring candidate directions and the same approach trajectory. The occurrence 
of double minima for this case may be noted. 

Figure 3. 2(j shows the fuel penalties for aim -point variations from the desired 

r . = 2838 Ion and i , = 116, 5°. The lower-right number at each intersection point 
pd d 

of the grid is the midcourse fuel required to target to that intersection, using the 

MFG law. The upper-left number at each intersection is the total correction fuel 

(midcourse plus trim) to trim the orbit to r and i . The values shown relate to 
h pci a 

a 2 midcourse correction to a slightly errant trajectoiy. The trim fuel values were 
determined with the numerical procedure for optimizing retro direction and firing 
time. It is apparent that the trajectory should be targeted to r ^ and at the mid- 
course correction. It may also be observed that the midcourse fuel cost for targeting 
to other grid points is very small. The fuel for trimming the trajectory targeted to 
r and i^ is seen to be about 3. 65 kg. It does not appear possible to separate trim 
costs due to inclination variations from trim costs due to radius of closest approach 

variations, since radius-correcting costs differ at different inclinations. Figure 3.2 7 

h 

differs from 3.2 5 only in that the midcourse maneuver is executed at 18 rather than 
at 2* 1 . Trim fuel costs differ somewhat between the two figures due to differing 
characteristics of the approach trajectory for the two maneuver times. Extra printout 
from one of the trajectories targeted to i = 106.5° revealed that the optimal firing true 
anomaly (for that trajectory) was about +16°, which would probably put the firing point 
closer to the equator than periapsis was. The retro direction and required trim fuel 
changed during the numerical optimization process as shown,. 



True Anomaly 

Right Ascension 

Declination 

Trim Velocity 

First Guess 

0 . 

82.81 

-41.77 

254.67 

Optimal Solution 

16.34 

117.19 

-48.41 

157.96 
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Figure 3.26 


Fuel Penalties for Aim-Point Variations 
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Analysis of in-plane retro and trim requirements is much simpler that for non- 
coplanar trims* Whatever the retro maneuver is, it usually results m an elliptical 
orbit about the Moon. If inclination errors are ignored, the optimal strategy for 
trimming to circular at a given radius is a two-impulse Hohmann transfer. Char- 
acteristics of such a transfer have been extensively studied. Figure 3.28 shows 
curves of trim velocity {sum of the two impulses) versus one apsis radius, called 
approach pericynthion. Each curve represents a different value of the opposite 
apsidal radius. The top curve is for 2338 km and the bottom curve is for 2838 km. 
Some of the middle curves are for opposite apsidal radii greater than 2838 km and 
some for smaller. Notice that the minimum of each curve occurs at no deviation 
from 2838 km in approach pericynthion and that the lowest curve also characterizes 
"correct” apsidal radius. Each of these cases requires only one impulse. Figure 
3. 29 was generated by contour-plotting the data of 3.28 . The x-axis of the plot 
is deviation of approach pericynthion from 2838 km and the y-axis is deviation of 
the opposite apsidal radius. The contour lines are separated by 10 m/s in trim 
velocity to circularize to 2838 km. The plot may be imagined to be an inverted 
pyramid. The value of such a plot in analysis of in-plane trim will be seen later. 


Trim and Retro Strategies 

Figure 3.3 0 is basically a contour plot of trim velocity in the plane of apsidal altitude 
errors. That is, the x-axis represents first-apsis altitude relative to 1100 km and 
the y-axis represents second-apsis altitude relative to 1100 km. The rectangles 
are lines of constant trim velocity to trim the orbit in-plane to 1100 km. They are 
spaced 10 m/s apart. 

Line is the locus of circular orbits. Points P^ represent possible pairs of first- 
apsis, second-apsis errors following a (hyperbolic) periapsis anti -velocity maneuver. 
Points located to the upper-left of L^ represent results of underburning with respect 
to circular at r and those to the lower-right of L^, overburning with respect to circular 
at r. The vectors emanating from P^ represent loci of points attainable by firing the 
retro at some radius greater than periapsis radius, always nulling the post-retro 
radial velocity with the retro maneuver. These loci are drawn as straight lines of 
slope 7. They closely approximate actual loci in the area of C = . 6 km 2 /sec 2 , 

O 
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Figure 3.23 


04/2O/T2 001 — i 

TRIM DELV VS APPROACH DISTANCE CPARM = OPAPSISD 



105 


Figure 3.29 
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r = 2838 km and dv = .7 km/sec. Let ns investigate each point. 
P 


P. and P. 
1 6 

p 2 

P andP, 

3 4 


The best strategy is to fire at periapsis. 

Fire at periapsis for best circularization but at the radius 
where the attainable locus crosses the y-axis for minimum trim. 

Better firing radii can be found for retro-fire for either circulari- 
zation or trimming to r^. 

The trim-minimum retro-fire location is periapsis, but the orbit 
could be circularized further out. 


* 


Final Midcourse Strateg ies 

Figure 3.31 shows the same trim contours as figure 3.30 , but has some different 

loci drawn over the contours. L„ and the two lines parallel to it are loci of points 

£ 

attainable for fixed C , 5v pairs by moving periapsis radius, r . It is reasonable to 

3 P 

assume that r can be changed by a midcourse correction with relatively little change 
P 

in C Q or dv. 


Any point in the half-plane above corresponds to an underburn at r^ while any point 
below L corresponds to an overburn at r . The figure shows that above L it is best 

Ci U u 

to make r^ = r^ so that the x-axis error is zero to minimize trim velocity (with no 
energy change). Below L , the midcourse maneuver should be targeted to the inter- 

Ci 

section of the r -locus and L . L is the locus of periapsis points (at various C *s 
P O o o 

and 6v's) for which it is possible to circularize with an in-plane retro maneuver at 
r^. This intersection is always to the left of x = 0. is the slope = 7 line passing 
through (0,0). If the last midcourse correction is targeted to r while the resulting 
energy and 6v correspond to a point below L 9 , the trim velocity will be minimized by 
retro-firing at a radius greater than r^ where y = 0 (second apsis is r^). 


3. 3.4 Pre-Targeting Results 

The pre-targeting capability of MAESTRO's midcourse guidance mode has been 

investigated. The results are quite satisfactory. Table 3. 7 is presented to 

illustrate the accuracy of the pre-targeting capability in predicting midcourse maneuvers. 

lx ix 

Minimum fuel guidance maneuvers were simulated in 10 steps beginning 10 after 
launch of a "hot" June 10, 1973 trajectory. The uncorrected trajectory is characterized by 
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The first line shown for each midcourse maneuver describes the pre-targeted (P) results. 

The second line represents the results obtained when the pre-targeted maneuver is 

applied impulsively to a multi -conic transfer. The third line shows the targeted results 

{still impulsive, multiconic). The constraint tolerances ar& 5 km in r and .2° in i. 

P 

The maneuver is described by impulsive velocity (m/s), right ascension (deg) and 

declination (deg). ARP is kilometers less than 2838 in closest approach radius. AT 

is degrees less than 116.5 in inclination. ATF is seconds less than 396000 (110 hours) 

2 2 

in time of flight. The energy, C^, is in km /sec . The first thing to notice from the 
table is how little the maneuver parameters and the flight time change during targeting. 


Table 3.7. 

Pre-Targeting Comparisons with Multi-Conic 
Minimum Fuel Guidance 


MC 

IT 

DVMG 

RAS 

DEC 

ARP 

AI 

ATF 

C3 

10 h 

P 

73.31 

72.88 

34.60 

- 

- 

6450. 

.630 


0 

73.31 

72.88 

34.60 

441.4 

2.70 

12670. 

.636 


1 

72.82 

73.02 

36.23 

.3 

.01 

12557. 

.637 

20 h 

P 

102,23 

77.85 

32.79 

_ 

_ 

10933. 

.643 


0 

102.23 

77.85 

32.79 

484.5 

4.85 

16455. 

.647 


1 

100.77 

79.46 

33.94 

1.8 

.03 

16931. 

.652 

30 h 

p 

129.92 

80.91 

31.87 

— 

_ 

14501. 

.656 


0 

129.92 

80.91 

31.87 

525.6 

6.72 

19217. 

.657 


1 

127.38 

82.35 

32.81 

2.3 

.04 

19665. 

.662 

40 h 

p 

160.13 

82.60 

31.34 



17279. 

.669 


0 

160.13 

82.60 

31.34 

506.1* 

6.65 

21041. 

.666 


1 

156.99 

83.47 

32.12 

1.7 

.03 

21206. 

.670 

50 h 

p 

197.35 

81.58 

30.91 

— 

— 

18317. 

.675 


0 

197.35 

81.58 

30.91 

493.1* 

6.63 

21087. 

.665 


1 

193.44 

83.12 

31.64 

3.3 

.05 

21635. 

.673 


* represents deviation from the desired radius of the variable approach 
guidance procedure 110 



Another item of interest in the table is the small size of the closest approach and 

inclination errors when the pre-targeted solution is propagated to closest approach 

with the multiconic method. The flight time for the pre-targeted (patched-conic) 

minimum fuel solution is seen to be about one or two hours longer than the targeted 

solution. The pre-targeted maneuver is compensated for this time difference by 

introducing a factor on the transfer time used in solution of Lambert's problem. 

A final item of interest in the table is the closeness of C g (Jacobian energy) for 

the pre-targeted solution and for the multiconic trajectory. The Jacobian energy 

is computed as the square of the Moon-relative velocity at patch to the sphere of 

2 2 

influence minus a constant term of about . 025 km /sec . The potential energy term, 

2 M (which is omitted from the energy computation in pre-targeting) has, by 
r . 

soi 2 2 

contrast, a value of about -. 148 km /sec . 


Tables 3.8 ■ and 3.9 are similar to 3.7 in that they show the progression of 
maneuver and end conditions during targeting. They describe FTA and VTE maneuvers 
at the same times as the MFG maneuvers of table 3.7 f and for the same errant 
trajectory. 


DVMG 

BTE 

BRE 

FTE 

VIE 

VCE 

TFE 

RPE 

INE 


midcourse velocity (m/s) 
error in B« T (km) 
error in B* R (km) 

flight time error (seconds less than 396000) 

hyperbolic excess velocity error (km/sec less than/762) 

circular excess velocity error (km/sec less than /jj/ ) 

r pd 


total fuel for in-plane trim to r^kg) 

radius of closest approach error (1cm less than r^ 2838) 

inclination error (deg less than 116.5) 


.3.5 Finite Burn Results 

Some comparisons have been made between the various methods for simulation of the 
midcourse burn phase. Of particular interest is how good an approximation the 
impulsive burn is to the carefully-integrated finite burn. Table 3. 10 . presents data 
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Table 3. 10 


Finite Burn Method Comparisons 


Method IT PVMG 


0 

0 

73.30 


1 

72.82 

1 

0 

73.30 


1 

72.90 

4 

0 

73.30 


1 

72.90 

6 

0 

73.30 


1 

73.03 


Method 

IT 

DVMG 

0 

0 

196.95 


1 

193.44 

1 

0 

196.95 


1 

195.97 

4 
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1 

195.66 

6 
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1 

195.04 
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436.4 

.4 

72.52 
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36.21 
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.4 

72.52 

73.15 

34.58 

36.21 

419.5 

.4 

72-. 52 
73.16 

34,58 

36.21 

400.5 

.3 
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31.64 

447.9 
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83.18 

30.98 
31. 63 . 

337.2 

2.4 
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for four methods of simulating the midcourse burn phase. 

Method 

Impulsive burn 
Cowell integration of the burn 
NIC E/Mean integration of the burn 
Closed-form/Multiconic burn approximation 

The trajectory represented in Table 3. 10 was a June 10, 1973 "hot” trajectory. 
Ignition times of 10^ and 50^ were considered. The 10^ correction required a burn 
duration of about 8.75 m while the 50* 1 correction required a burn of about 31. 5 m 
(too long, considering available fuel) for the minimum fuel maneuver. The zero-th 
iteration for each method was primed with the same impulsive from pre-targeting. 
The pre-burn and post-burn trajectory method was multiconic in each case. The 
column headings have the following meanings. 

IT Iteration number 

DVMG Impulsive velocity magnitude (m/s) 

RAS Right ascension of thrust (deg) 

DEC Declination of thrust (deg) 

RPE Radius of closest approach error (km less than 2838) 

INE Inclination error (deg less than 116.5) 

FTE Flight time error (seconds less than 396000) 

VIE Hyperbolic excess velocity error (km/sec less than /T62) 

TBURN Midcourse burn duration (minutes) 

* 

Each midcourse burn is targeted to RPE and INE tolerances of 5 km and .2 deg, 
respectively while constraining the maneuver AV to lie in the critical plane. Several 
observations can be made about the data presented in Table 3.10., 

1. The impulsive approximation errors are only about the size of expected 
resolution and pointing errors for the earlier correction and are still 
almost negligible for the later correction. 


IBURN 

0 

1 

4 

6 
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2. the variation of the Method 6 solution from the impulsive solution 
is greater than that of Methods 1 and 4 for the 10 h correction, but 

i 

less for the 50 correction. 

3. Flight time for the MFG solution does not change very much between 
methods . The change in FTE between impulsive and other methods 
is much less than the burn duration, indicating suitability of applying 
the impulse at ignition time. 

L. 

4. Methods 1 and 4 give exactly the same results for the 10 correction, 
but differ for the 50 h correction. The difference may be due to the 
integration step size being too large (it is the burn duration) for one 
method or the other. 
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SECTION 4 

ON-UNE CONTROL AND DISPLAYS 

This section will describe the methods for communication between the Flight 
Director and the mission control program. The Flight Director must have a means 
of controling the scope of the analysis. He will determine which analysis will be 
performed at which times. The Flight Director must also be provided with displays 
which present the information in a concise, easy-to-read format. The first section 
describes the method which will enable the Flight Director to control the analysis. 

The second section describes the displays which are generated for the Flight Director. 

4.1 CONTROL 

Figure j s an overall block diagram of the mission control system. The system 

is controlled by the Flight Director. He is in charge of all mission analysis operations 
and is responsible only to the Mission Director. The Flight Director will control the 
analyses through the 2250 console using the Graphic Terminal System. GTS is a system 
executive which controls the submittal of jobs to the 360/75 computer. This system 
also has the capability to retrieve and edit data files with the use of the 2250 light pen 
and typewriter. By using this system, inputs which determine the type of analysis 
performed by MAESTRO can easily be controlled by the Flight Director. GTS has been 
used during several simulations and has functioned satisfactorily. The use of GTS will 
not degrade the reliability of MAESTRO because the card input capability of MAESTRO 
is still retained. Thus, MAESTRO can still be used in the event of a failure of the 
2250 or on another computer if scheduling conflicts exist with the orbit determination 
or attitude determination programs on the 360/75. 

The Flight Director has the responsibility for all operations of MAESTRO and he will 
decide which analyses are required at which times. A team of spacecraft experts and 
the mission scientist are available for consultation with the Flight Director. This team 
will help evaluate the effects of mission analysis decisions on the basic scientific 
objectives of the mission and operational characteristics of the spacecraft. 
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Mission Control System Block Diagram 



Figure 4.1 



























MAESTRO was designed with many operational modes to provide the Flight Director 
with the capability of performing the necessary mission analysis tasks as rapidly as 
possible. The operational modes of MAESTRO were defined in Section 2. Here they 
will be described according to their roles in the decision-making process. 

A. Preliminary GO/NOGO analysis 

As soon as a preliminary state is available from the orbit determination 
program, midcourse analysis runs will be made to determine the amount 
of fuel required and approximate execution time for the first MC maneuver. 
This data provides the Flight Director with the information to determine 
whether the primary mission is achievable and the amount of analysis time 
available until the first command decision must be made. 

B. Refined Midcourse Analysis 

Many midcourse analysis runs will be made with different guidance laws 
and execution times to determine the final midcourse velocity and attitude. 

C. Verification of the Midcourse Decision 

Precise midcourse and approach analysis runs will be made to verify the 
decision made in B, Monte Carlo runs will be made to verify that the 
decision will provide a high probability of mission success. 

D. Lunar Approach Analysis 

Approach analysis runs will be made to determine the retro firing time 
and attitude. 

E. Lunar Approach Verification 

Precise approach analysis runs will be made to verify the decision made 
in D. Monte Carlo runs will be made to verify that the decision will result 
in a high probability of mission success. 

F. Post Injection Trim 

The post injection trim analysis will be used to determine the firing times 
and attitude of the impulse(s). This analysis will be used in both an approx- 
imate and precise modes. 
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G. Alternate Mission 

In the event of a launch vehicle of spacecraft failure it may be possible to 
achieve an orbit which will still yield valuable scientific data. The midcourse 
or approach analysis modes can be used to determine the maneuver to achieve 
an acceptable alternate mission. These analyses can be performed in both 
approximate and precise modes. 

4. 2 DISPLAYS 

MAESTRO will have the capability of presenting displays via the 2250 console, on 
hardcopy or as the output from a plotting processor. These displays will have levels 
of complexity that will range from an easy to read outline form to a comprehensive 
presentation of the analysis- performed. The current displays available in MAESTRO 
are shown in the uses guide, reference 2 . These displays are obtained on the 2250 
console under GTS and hardcopies are also printed. It was found that the 2250 screen 
is hard to read because the images were constantly blinking. Thus, hardcopies of the 
2250 displays should be provided to the Flight Director as soon as possible. 

The displays contained in MAESTRO provide the Flight Director with a means for 
rapidly accessing the spacecraft's capability under various in-flight mission situations. 
The basic ground z’ule used in designing the MAESTRO displays was to provide the 
Flight Director with all of the information he needs without clouding the problem with 
too much information. In some cases this necessitated more than one display for an 
analysis mode. For example, the approach analysis displays have two levels of com- 
plexity. One display provides the Flight Director with a rapid estimate of the retro 
motor's attitude which results in a minimum post injection trim requirement. A second 
display provides comprehensive information about each trial firing at the desired 
attitudes. All of the displays available in MAESTRO are discussed in the following 
paragraphs. 

RETRO FIRE ANALYSIS 

There are two different displays presented in this analysis. One display is a summary 
type of display presenting minimum eccentricity and post injection trim fuel requirement 
at the retro attitudes tried. This display is actually a contour map printed by the 

computer and depicts the attitudes which result in the minimum trim fuel requirement. 
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The Flight Director can use this display to isolate the desired range of attitude. 

A second display presents, detailed information about each trial retro firing. This 
display contains the mission constraints and post retro orbital elements as a function 
of firing time for each attitude. The display is set up so that a page is devoted to 
each retro attitude. Thiis, it is easy to see that the Flight Director could easily be 
snowed under by the output from this mode. The first display is included to eliminate 
the undesirable attitudes so that the Flight Director will not have to waste Ms time 
sifting through a large stack of output searching for good attitude angles. 

The information from the retro firings can also be saved on a tape. This tape can 
later be used with some type of post processor or a plotting program, 

MID COU RSE ANALYSIS 

The displays presented by the midcourse analysis provide the Flight Director a means 
of accessing the spacecraft's capabilities as a function of the available control variables. 
Logic is included to automatically scan the execution time and flight time so that a 
display presenting the fuel requirement and constraints as the above quantities are 
varied can be presented. The display also contains information about execution errors 
and retro conditions. The execution error portion of the display contains the error 
elipse at the target planet of the quantities targeted. It also contains the expected 
second midcourse requirement. The retro portion of the display presents the retro 
firing time, attitude, mission constraints and post injection trim requirement after 
retro for each midcourse correction. The quantities derived for this display assume 
an errorless midcourse correction, 

A second display is also available from the midcourse analysis mode. This display 
presents the target planet conditions as a function of firing time, midcourse velocity 
and attitude, TMs display was included for two reasons. First, as a failure mode 
analysis tool in case of attitude control failure. The midcourse requirement can be 
determined for any fixed attitude or range of attitudes if the spacecraft's attitude is 
constrained for some reason, 

TMs display can also be used to determine the stability of some predetermined mid- 
course maneuver. The execution time, velocity and attitude can be scanned about 
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some selected maneuver to determine the affects of errors of one of the above quantities 
on the conditions at the target planet. 

A tape is also written from this analysis for use in a plotting program. 

MONTE CARLO ANALYSIS 

The Monte Carlo analysis display contains information about each sample and a statis- 
tical summary at the end of the analysis. The information printed for each sample is 
the fuel, its mean value of the past samples for the first and second midcourse, first 
second and third trim maneuver and fuel remaining. The spin axis-sun angle, its mean 
and standard deviation is also presented for the maneuvers mentioned above along with 
the retro spin axis sun angle. The pre-trim orbital elements, their mean and standard 
deviations are also presented. These elements include the apo and periapsis radius, 
argument of perigee, inclination, longitude of the ascending node, true anamoly and 
retro firing time. The information presented for each sample provides the Flight 
Director with the specific details of each sample and indicates when enough samples 
are run. As an option, the Flight Director can also obtain information to restart 
any sample which he things might require more analysis. 

A statistical summary is printed at the end of the display. This summary presents 
the mean, maximum value encountered and standard deviation of all the quantities men- 
tioned above. This includes the fuel and spin axis-sun angle of all the maneuvers and 
the pre-trim orbital elements. 

VERIFICATION ANALYSIS 

The display contained in the verification analysis can be divided into three parts, real- 
time thrusting analysis, station coverage and target planet conditions. The real-time 
thrusting display is used to obtain information that can be used while the midcourse 
motor is burning. This display presents the spacecraft thrust, weight and velocities 
away from the visible tracking stations as a function of time since ignition assuming 
a nominal engine burn. The thrust data can be converted to chamber pressure data 
and an expected doppler shift determined from the velocities away from the visible 
tracking stations. This data can then be compared to the real-time telemetry data 
and corrective action initiated if the telemetry data wildly differed from the predicted 



The tracking coverage display presents the local time* elevation angle and azimuth 
angle of the spacecraft for each of the tracking sites at the time of midcourse engine 
ignition. 

The last part of the verification display presents a summary of the motor firing con- 
ditions and the conditions at the target planet. The motor conditions include the burn 
time, attitude and spin axis -sun angle at ignition. The conditions at closest approach 
to the target planet are presented. These are the time of closest approach, radius 
and inclination with respect to the target planet’s equator, 

POST I NJ ECTION TRIM ANALYSIS 

The post injection trim display available presents the elements of the initial orbit, 
transfer orbit and final orbit. The firing time, magnitude, direction and spin axis- 
sun angle of each trim maneuver is also output. This display is automatically presented 
for two types of trim maneuvers. First, a two impulse plane-change maneuver. 

The other type of maneuver is a two impulse Hohmann transfer maneuver with no plane 
change. 
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SECTION 5 

OPERATIONAL LOGISTICS AND TESTING 


This section contains a description of important factors concerning the decision- 
making processes. The first subsection presents a typical sequence of events . 
for this mission and emphasizes the necessity for a smooth man-machine 
interaction. The second section describes possible failure modes and presents 
contigency plans in case some of these failures should occur. The last section 
shows how tests have been conducted on the MAESTRO system to insure its 
reliability. * . 

5. 1 FLIGHT PROFILE OUTLINE 

In Reference(l) a fairly detailed outline of the probable sequence of events was 
presented but there was little said about the timing of events. Some preliminary 
estimates on the timing for the key elements of the system are now available 
and it is important to consider the sequence in which these elements will be 
used. 

Suppose, for clarity, that the launch is initiated at Cape Kennedy on June 10, 1973, 
12:00 hours GMT, Suppose further that the nominal third stage burnout time is 
15 minutes after liftoff. 

At the time of liftoff, or before, the entire mission control system should begin 
a test exercise to verify that all elements of the system are operational at that 
time. All data channels, storage devices, and operational modes of the system 
should be briefly tested in a way that will show up any possible problem. 

Let us now suppose that the third stage injects the spacecraft into its translunar 
trajectory at 12:15 GMT and that tracking data begins coming in a 12:30. This 
possible delay is trajectory dependent but w r e assume for illustration that the 
time between translunar injection (TU) and first acquisition by one of the 
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trfinkircr stations is 15 minutes. Assume that it requires 15 more minutes of 

" ■* • • o 

tracking, data processing, and differential correction to obtain a reasonably 
accurate but unrefined state vector. This vector may now be entered into the 
mission control program for a preliminary midcourse sweep 1) to establish 
a rough idea of the situation and 2) to dtermine the approximate amount of time 
available for analysis and attitude maneuvers before the first midcourse 
correction. 

At the same time tracking data is first acquired, telemetry from the attitude 
sensors can also be expected. The relative speed with which the orbit and 
attitude can be determined may be more important than assumed here but let 
us suppose that a rough estimate of the burnout attitude can be obtained from 
the direction of velocity vector at TLI or from the nominal value of that 
direction. In any case, it is time for MAESTRO to determine the time and 
magnitude of the required 1st midcourse correction. 

The midcourse sweep mode of the MAESTRO can be run in about 30 seconds 
(360/75) per guidance law with 10 trial firings ranging from 2 to 20 hours past 
liftoff. After the initial rough determination of the state, the major resources 
of the computing system should be turned over to the MAESTRO while three or 
four guidance laws and an approximate approach analysis are run. The time 
required will be about five minutes and the computer should then be instantly 
returned to the Orbit and Attitude Determination Programs (ODP and ADP) for 
intensive tracking and the refinement of the estimated state. 

The early stages of control, then, can be summarized as follows: 

GMT 

12:00 Lift off-start of final testing operations 

12:15 Launch phase completed 

12:30 Acquisition by Tracking System 

Initiation of Orbit and Attitude Determination Operations 
12:45 Preliminary Anchor Vector Ready 

Possible Rough Attitude Data 
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12:50 


Preliminary Midcourse Corrections and Approach Analysis 
Initiation of Intensive Orbit and Attitude Determination 
operations 

Probable Acquisition by 2nd Tracking Station 
Evaluation of Midcourse Sweep Results 

Assume, now, that the energy of the translunar orbit is considerably in error 
and that the midcourse sweep results indicate that the 1 st midcourse correction 
should be performed as soon as possible. The attitude maneuver sequence 
should be initiated immediately so as to permit the early correction.- Suppose 
that the decision to begin the maneuver is taken in the 10 minutes between 12:50 
and 13:00. The remainder of the control sequence might be as follows: 

13:00 Accurate Attitude Available 

At least two stations tracking spacecraft 
Decision to initiate attitude maneuver 
13:05 Commands to Spacecraft Attitude Control System 

ODP has computer for refined orbit determination 
13:15 Refined orbit available 

Attitude maneuver in progress 

Computer to MAESTRO for detailed midcourse analysis 
13:20 Attitude maneuver in final stages 

Detailed midcourse analysis complete 
(De sired Attitude Known) 

Final trim pulse of ACS 

Computer to ODP and ADP for verification 

The haste with which the above operations have been performed correctly 

indicate that the hypothetical situation is critical and that further delay will 
probably result in an unsuccessful mission. At this point, then, a critical 
decision is required of the mission and flight directors. If the ODP and ADP 
have relatively good estimates of the state, it may be wise to dispense with any 
refinements in the control process and to perform the first midcourse maneuver. 
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The decision-snaking process would be a trying one in which the risks of making 
a hasty decision would have to be carefully weighed against the probability of 
having insufficient fuel to complete the mission. 

The consideration of such a rapid 1st maneuver control sequence assumes that 
very unlikely launch vehicle dispersions have occurred but it seems important 
to determine the worst possible situation that could still be turned into a 
successful mission by emergency use of the midcourse correction system. 
Further, there may be good reasons, of which this writer is unaware, why 
the initial control sequence outlined above simply cannot be executed- in the 
imaginary time scale of the outline . The outline, then, should be considered 
as a suggestion that the earliest time at which an emergency 1st midcourse 
maneuver can be initiated is approximately 90 minutes after liftoff. The 
schedule above is, in any case, a goal to aim for in establishing the response 
time of the system. In most cases, the much more conservative outline of 
the design study final report should be followed. 

We come now, to the consideration of the maneuvers to be carried out after 
the spacecraft has been placed in lunar orbit by the solid retro-rocket. 

A very important difference in the situation will now be realized; the operations 
will no longer be time-critical. In the event of a large retro-rocket errors or 
other complications, the mission may be fuel-critical but we can expect to 
have plenty of time for analysis and should not unduly rush the trim operations. 
Indeed, it will be desirable to track the spacecraft for some time after lunar 
orbit insertion in order to study the effects of the lunar gravity on the orbit 
evolution and, thereby, to verify the adequacy of the existing models of the 
lunar potential. 

Consider then, as we always should, what might happen in the real-time 
control center after the news is received that the spacecraft has been successfully 
placed in lunar orbit. The optimal orbit insertion attitude (from dynamic 
considerations) will probably not be compatible with the desired cruise-mode 
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spin-axis orientation and it will probably be necessary to reorient to 
an attitude that is commensurate with the power and communication 
requirements of the spacecraft. Of course, the lunar orbit should be 
checked for short-term stability before a reorientation is performed. 

This check will guard against an early impact or wild gyrations of the orbit 
in case of very large deviations from the nominal retro performance. The 
check, however, should be an integral part of the lunar orbit insertion 
verification procedure and is not considered as part of the time-uncritical 
trim maneuvers. 

The post-insertion operations, then, should begin with an orbit verification 
procedure and the operations should approximately follow the outline 
below. 

I. Lunar Orbit Verification Procedure (Time Critical) 

a. Post Injection Orbit Determination. 

b. Post Injection Attitude Determination. 

c. Verification of Short-Term (several weeks) Orbit Stability 

1. Stable Orbit - Go to n 

2. Unstable Orbit - Commence Emergency Trim Maneuvers. 

II. Attitude Reorientation Decision (Time Uncritical) 

a. Existing Attitude Suitable for Cruise Mode - Go to HI. 

b. Existing Attitude Unsuitable for Cruise Mode - Reorient 
to Desired Cruise Mode Attitude. 

HI. Verification of Lunar Gravity Models 

a. Intensive Long-Term Tracking. 
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b. Existing Gravity Models Adequate for Long-rerm Orbit 
Prediction - Go to IV. 

c. Existing Models Inadequate - Perform Lunar Gravity Model 
Determination {separate program - not in Mission Control program) 

IV. Trim Maneuver Policy Decision 

a. Number of Maneuvers Required 

(Minimum Fuel vs. Risk of Multiple Maneuvers) 

> 

b. Effects of Shadowing and Occultation on Trim Policy. 

V. Statistical Analysis 

a. Probability of Mission Success vs. Trim Motor Ignition 
Time and Attitude. 

b. Definitive Orbit Evolution Prediction. 

VI. First Trim Maneuver 

a. Reorient Spacecraft for Trim Motor Ignition. 

b. Commence Trim Motor Ignition Sequence: 

c. Monitor Telemetry During Trim Maneuver. 

VII. Post-Trim Evaluation 

a. Intensive Calculations. 

b. Orbit Evolution Calculations. 

c. Mission Requirements Satisfied? 

No - Go to IV and repeat. 

Yes - Turn over spacecraft to boom - deployment experts. 
Flight dynamics control completed. 
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5.2 


CONTINGENCY ANALYSIS 


It should be noted at the beginning of this section that the items of discussion 
are necessarily of a pessimistic nature and that there will probably be a 
tendency for the people involved to rebel against the consideration of the 
many possible ways in which the mission might fail of partially fail. No 
one wants the mission to fail and the discussion of the possibilities may set 
loose psychological factors that will be detrimental to the objective consideration 
of the action to be taken in such situations. 

For example, if we consider the possibility that the computer logic in the 
mission control program will fail in some unforseen situation, the planners 
and programmers may have a tendency to discount the probability of occur- 
rence of such difficulties. If we discuss the possibilities of the failure or 
partial failure of the spacecraft hardware, the designers of the spacecraft 
system may feel slighted and unwilling to consider some failure modes as 
viable possibilities. It seems important to attempt to dispel such inclinations 
at the outset so that effective contingency measures may be preplanned in a 
way that will not be clouded by the egocentric factors that affect us all. 

Let it, therefore, be clearly understood that the following suggestions of 
possible failure are made with purely constructive intent. 

The major areas of possible failures can be classified as follows: 

I Hardware Malfunctions 

n Computer Software Errors 

HI Human Errors 

Although the Hardware and Computer Software categories might, in a broad 
sense, be considered human errors, the third category now refers to those 
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errors in judgement or communication during the in-flight operations that 
could, in post-flight analysis, be considered causes of an unsuccessful 
mission. 

I Hardware Malfunctions 

This category is the most extensive because its possibilities affect the entire 
in-flight operation to a marked degree. The malfunction of any piece of 
hardware opens up many possibilities for computer software inadequacies 
and human error, where- as in-flight failures in categories n and in do not 
affect the operation of the spacecraft. 

The possible hardware malfunctions are numerous and are further class- 
ified as follows: 

A. Spacecraft Hardware Malfunctions 

1. Launch Vehicle (large Errors) 

2. Midcourse Correction and Lunar Orbit Trim System 

3. Attitude Determination and Control System 

4. Lunar Orbit Injection Retro Rocket 

5. Spacecraft Power and Communication Systems 

B. Computer Hardware Malfunctions 

1. Input/Output Hardware Failure 

2 . Partial Memory Failure 

3 . Mainframe down temporarily 

4. Power Failure 

C. Tracking System Malfunctions 

1. Communication link failure 

2. Station malfunctions 

3. Data Processing Equipment failure. 

Di Uplink Command Hardware Malfunctions 
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Computer Software Errors 

A. Orbit Determination Software 

1. Preprocessing Software Failure 

2. Diffential Correction Program Failure 

3. Data Transmission or Storage Errors. 

B. Attitude Determination Software 

1. Preprocessing Software Failure 

2. Attitude Calculation Errors 

3. Underestimate of Required Attitude Control Fuel / 

4. Data Transmission or Storage Errors 

C„ Mission Control Program Software 

1. Midcourse Guidance Calculation Errors 

2. Required Attitude Calculation Errors 

3. Retro Rocket Ignition Time Calculation Errors 

4. Incorrect Evaluation of Required Trim Fuel 

5. Data Transmission or Storage Errors 

Human Errors 

A. Errors in Planning 

1. Errors in assignment of priorities during the flight 

2. Errors in scheduling work 

B. Errors in Communication 

1. Incorrect or Poor Descriptions of Situation 

2. Unclear Definition of Work Assignments 

3. Unclear Presentation of Results 

4. Poor Definition of Objectives and Priorities 

C. Errors in Judgement 

1. Incorrect Evaluation of the Situation 

2. Faulty Evaluation of Probability of Success 



The outline above is a compilation of the areas in which malfunctions and 
errors might occur that would be either totally or partially disastrous. 

There may be other important possibilities that have not been considered. 

In the discussion that follows, three of what are thought the most likely 
failure situations are considered. 

Three important failure modes that could cause the mission to be unsuccessful 
are: 

a. Midcourse and Trim Propulsion System Failure 

b. Attitude Control System Failure 

c. Retro Rocket Failure 

These three most drastic possible malfunctions are considered below in 
order to suggest what might be salvaged in the event of such failure, 

a. Midcourse and Trim Propulsion System Failure 

The malfunction or failure of the midcourse system at any time 
prior to lunar orbit insertion would turn the mission control problem 
into an AIMP-E type decision-making process. The attitude and 
ignition time for the retro-rocket would then constitute the only 
control parameters available for use. 

The current capabilities of the MAESTRO are adequate to provide 
sufficient information for the mission director to determine the "best" 
attitude and retro ignition time from the standpoint of the mission 
requirements. The approach analysis and retro attitude optimization 
are integral parts of the current system and could be used in the 
deterministic or Monte Carlo modes of the program to aid the mission 
director in making the probably unsatisfactory "best" decision in case 
of midcourse and trim system failure. 

b. Attitude Control System Failure 

This possible problem has been considered and steps have been taken to 
permit inflight evaluation of such a contingency. The solution, in the case 
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of ACS failure, would be to determine the ignition times and 
durations of the many midcourse corrections that would be 
required to salvage the mission. 

There seems to be no reason why, with the large propulsive 
capability of the spacecraft, a satisfactory lunar orbit could not be 
achieved in many cases. Of course, the results would depend upon 
when the ACS failed (that is, upon the actual value of the fixed attitude) 
but it seems obvious that many possible situations could be handled 
in a satisfactory way be multiple midcourse corrections and trim 
maneuvers . The capability to analyze such situations has been 
incorporated into MAESTRO in the form of a fixed attitude guidance law. 

c. Retro Rocket Failure 

This possibility is probably the most catestrophic of all because of 
the large energy change supplied by the (single) firing of the solid 
rocket. 

It is this writer's understanding that some of the experiments on-board 
the RAE-B spacecraft would not be seriously impaired if the orbit were 
highly elliptic with respect to the Earth. These experiments might be 
saved or helped by the use of the midcourse correction system after 
the retro rocket failure to prevent a very close approach to or impact 
on the Earth. 

When the flight plan called for a leading approach to the Moon, it was 
conceivable that the midcourse and trim propulsion system could be 
used, after jettison of the faulty retro rocket, to adjust the highly 
elliptic Earth orbit so that a substantial lifetime could be achieved. 

Since the revision of the flight plan to yield a trailing approach, it has 
been found that the nominal trailing swingby precludes the possibility of 
using the midcourse motor to establish a stable Earth orbit. The 
. energy gained during the close trailing approach is so great that the 
spacecraft will escape the Earth-Moon system unless the retro- 
rocket is fired and the midcourse propulsion system cannot provide 
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enough impulse to reestablish an elliptic orbit about the Earth. 

The failure of the retro-rocket, then, should be considered 
catastrophic for the mission if the lunar approach is near 
nominal. If, however, there is some early indication that 
the retro-rocket might not fire, it may be wise to change the 
aiming point so as to lower the energy gained during the flyby so 
that the spacecraft will not escape the Earth-Moon system. Studies 
are underway to evaluate the cost and the Earth orbit options 
available in case of this not very likely situation. 

> 

5.3 SYSTEM TESTING 

MAESTRO was tested in December, 1971 during two three-hour sessions on the 
GSFC 360/75 computer. The program was operated under control of the Graphic 
Terminal System via the 2250 graphic display console. The various options of 
the program were tested for a simulated in-flight situation. Unfortunately, 
due to scheduling problems and computer hardware failures, it was not 
possible to exercise the program interface with the Orbit Determination 
Program but some tests were made in which the ODP and MAESTRO were in 
core at the same time. It was found that either of these highly compute bound 
programs can cause the other to run considerably slower than it would otherwise. 
If the 360/75 is to be the primary computer for the launch, it will be 
necessary to arrange a rigid schedule for assignment of priority to the two 
programs. In the event of a large translunar injection error it may be necessary 
to rim the two programs on separate computers in order to insure that the first 
midcourse calculations can be performed in time to make the actual maneuver. 

The internal operation of MAESTRO during this simulation was quite satisfactory 
although several minor problems were discovered and corrected. Had the 
spacecraft actually been in flight, it would have been possible to execute the 
first midcourse in a way which, although not optimal, would have resulted in 
a successful mission. Assuming that the first midcourse was executed correctly, 
it would have been possible to determine the retro firing time and 
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attitude which would result in a successful lunar orbit. 

The long-delayed first full -participation simulation was held on Wednesday, 

June 14 and the results were very encouraging. With a few noticeable exceptions, 
the simulated decision-making process went very smoothly and provided a valuable 
first step in the consolidation of the three major components of the system. It was 
obvious that everyone had done considerable homework in preparation for the 
exercise and that most of the problems were do to interactive aspects of operation 
father than to lack of preparation. The remainder of this section will be concentrated 
on the problems so as to focus attention upon them and prevent their future occurrence. 


The major operational difficulties can be classified in three areas; 

1. System Problems (GTS) 

2. Communication Problems 

3. Programming Problems 

The following is an abbreviated record of the simulation in chronological order: 


21.00 GMT 
21.11 


21.24 

21,37 

21.41 


June 14, 1972 
360/75 available for use 
Graphic Terminal System (GTS) up 
Two Problems 

1. Wrong version of GTS 

2. Attitude Determination Program (ADP) not ready for 
loading 

Attempted Reloading of GTS and programs 
Unsuccessful - JCL error 
All programs ready for operation 

Tracking data (simulated) coming in - GPUT processing data 


21.43 Orbit Determination Program (ODP) processing data 

21,48 Two iterations of differential correction (DC) process 

complete in ODP 

ADP processing (simulated) data using nominal orbit data 

21, 55 Flight Director (FD) waived 4th iteration of DC. ODP out, 

MAESTRO(M) in core. 
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22.02 Attitude available to within about 3 degrees assuming 

nominal orbit. 

ADP was not aware that M required attitude information 
at this time, 

22.04 M starting midcourse guidance calculations with anchor vector 
(AV) from 3rd iteration of DC. 

22.08 Fixed Time of Arrival (FTA) Guidance sweep completed 

Problems: 1. Midcourse guidance calculations taking 
4 times longer than expected. 

2. Rewind error in M - corrected through input array. 

22.12 M out - ODP in core 

GPUT processing second pass of simulated data. 

22.19 Problems: 

1. ADP unable to read AV from ODP 

2. ADP previously unaware that 1st AV was available 

AV entered into ADP by hand. 

22.20 ADP running with 1st AV 

22.25 Attitude available to about 1° 

22.26 ODP beginning to process second pass of simulated data 

22.42 First iteration of 2nd pass data finished on DC 

22.43 Second AV to M - ODP out of core. 

Start midcourse sweep with 4 guidance laws. 

22.46 ADP unaware of new AV 

22.57 Midcourse sweeps with 4 guidance laws completed on M. 

Problem: Slow calculations confirmed. Midcourse sweeps 
required 12 minutes - earlier tests with same software 
showed 4 minute running time. 

22.58 New AV to ADP 

23. 05 M still in core - Beginning Midcourse Verification at 10 hours 
with FTA guidance law. 

23.10 Midcourse verification still running - should have finished 
in a few seconds. 

23.11 Midcourse run with numerical integration submitted. 
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23.12 


ODP operators confirm very slow running time problem 
ratio 4 or 5 to 1. 

23.21 Midcourse verification run lost. Very large I/O time for 

run in accounting information. 

23.27 Desired attitude to ADP for Attitude Control System (ACS) 

command sequence calculations. 

23.30 M out, ODP in core for simulation of lunar orbit tracking. 

23.32 Bad start of ODP - incorrect input caused termination. 

23.33 ODP operational again. 

23.41 ADP finished with first try at generation of ACS command 

sequence. 

Problem: Command sequence incorrect. 

Timing for command sequence approximately correct. 

23.47 Two passes of lunar orbit data to ODP - starting DC, 

24. 04 Starting 7th iteration on DC for lunar orbit - results fair- 

trying one more iteration. 

24. 10 Simulation halted. 

The outline above shows the general flow of the simulation and the problems that 
were encountered. Of the most importance are the timing of the separate operations 
and the interaction (or lack of it) of the major components of the computing system. 

After this successful June simulation, testing efforts were devoted mainly to 
preparation for the IMP-H launch operations in September 1972. These tests 
were conducted with a modified version of the MAESTRO program and , while 
the attitude and orbit determination programs were different from the ones 
that will be used for RAE-B, the required operations were very similar to those 
for the lunar orbit mission. The use of MAESTRO for IMP-H was, by far, the 
best kind of testing we could have wished. The differences, both psychological 
and operational, between a simulation and an actual launch are so great that the 
experience gained during the IMP-H launch was worth several RAE-B simulations 
in spite of the differences in the flight prof iles. 
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Since the September launch of IMP-H, there have been several attempts to conduct 
a realistic RAE-B simulation but these have included little more complexity than the 
testing of the interfaces between the various computing systems. Very recently, some 
more thorough simulations have been held but the author was not involved and cannot 
report on the effectiveness of the tests. A significant feature of the most recent 
tests is that simulated data can now be generated during the simulation and this 
capability adds a great deal of flexibility to the operational tests. 

A schedule has been established in which simulations are to be held every other week 
from now until the launch. Every effort should be made to make these simulations 
as realistic as possible and a standardized schedule of operations should be established 
as nearly as is possible with the changing computer hardware and software. 

There is an important aspect of the entire MAESTRO design and development work 
that has received no overt recognition but has contributed to the effectiveness of the 
preparation for the launch. This aspect is the freedom to criticize the launch support 
operations as the shortcomings become noticeable. It is not usually acceptable for 
a contractor to criticize the Government but, throughout the MAESTRO development, 
there have been several points which have been the subject of controversy and for 
which our opinions have been freely voiced with the most constructive intent. This 
freedom to discuss what is wrong has been very helpful in allowing us to come to 
grips with many problems and has increased the usefulness and launch readiness 
of the MAESTRO program. 

Internal checks are constantly being conducted with MAESTRO. This is accomplished 
by comparing the results of the various program modes against each other and 
against independent programs. The use of the GSFC 360/91 computer for program 
development has greatly assisted in the testing of MAESTRO. Goddard personnel 
are now able to use new program capabilities as soon as they are available or, 
sometimes, while still in the development stage. The error detection abilities 
of Goddard's personnel are uncanny. Their help in this area and their general advice 
is greatly appreciated. 
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SECTION 6 


CONCLUSIONS 


A description of the MAESTRO program has been presented. The capabilities 
and various operational modes have been described in a general way and the 
reader has been referred to the programmer’s and user’s manuals for more 
details. It has been emphasized that the program's flexibility is one of its 
principal assets and that the all-FORTRAN block-style development is largely 
responsible for the internal flexibility and ease of modification. 

The major operational* modes are: 

1. Trajectory Propagation 
Method and Integrator Optional 

2. Midcourse Guidance or Kick-Motor Firing Analysis 

Includes Statistics of Propagated Errors 

3. Lunar or Planetary Orbit Insertion Analysis 

4. Orbit Trim Analysis 

One, Two, or Three Impulse Transfer 

5. Midcourse and Orbit Trim Verification Analyses 

•Simulations of Actual Firings 

6. Orbit Evolution Analysis 

Uses numerical averaging or methods of (1) above 

7. Monte Carlo Analysis 

Uses any of above in realistic mission sequence- 

In addition to these capabilities, a modified version of the program has been 
developed that has the capability to perform parameter estimation analysis 
to recover the gravitational harmonic coefficients for the central planets' field 
as well as parameters that describe solar radiation pressure effects. This 
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parameter estimation version of the program features a numerical integration 
of the linear variational equations for the method 7 elements or an integration 
of the averaged linear variational equations a nalogous to the method 8 
numerical averaging. The parameter estimation version is described in 
Reference (5) which supplements this report. 

It has been pointed out that MAESTRO was designed to be a multi-purpose 
mission analysis tool as well as an in-flight mission control program. The 
program is capable of use for design of Earth, Moon or planetary orbit or 
flyby missions and is capable of supporting the mission control operations 
for such missions. The logic for display of mission-peculiar information is 
isolated from the program options and it is a relatively simple matter to 
modify the output for any specific spaceflight mission. It is recommended 
that care be taken to preserve this isolation and flexibility if the program is 
modified for use on other missions. 

Some care has been taken to make the program as machine independent as 
was commensurate with the required interface logic for direct communication 
with the attitude and orbit determination programs. It is suggested that any 
future modifications be included so as to maintain this flexibility. It is hoped 
that the MAESTRO program and the techniques described in this report will 
be useful in the future exploration of space. 


140 



APPENDIX A 


INTEGRATION OF MODIFIED ELEMENTS 


The terms to be integrated are, 

1. semilatus rectum, p 

2. e sin to 

3 . e cos to 

4. to+ f 

5. inclination, i 

6. longitude of the ascending node, D 

where e is the eccentricity 

w < is the argument of perigee 
f is the true anomaly. 


The derivatives of the above quantities are numerically integrated to obtain the 
instantaneous orbital elements. The derivatives of these quantities are derived 
in many texts and reports and only the results are presented. The derivatives of 
the elements are 


= I 2 r 9 


( e sin to) 



■V? 


- cos (to+ f) R + 
r 


[M - 


sin (io+ f) + —~ e sin to 

Jr 


- e cos to — sin ( to + f) cot i W 


(e cos to) = 


JZ- 


sin 


(to + f) R + J^l + cos (to+ f) + e cos toj 3 


+ — e sin to sin (tu+ f) cot i W 
P . 


(to+ f) = /r 2 - sin (^+ f) cot ij 


W 


1 = 

« 


r cos (w + f) W / y p r u 
r sin ((,•+ f) W / (sin i 
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where ^ is the gravitational potential 
r is the radius 

R, S, W are the components of the inertial perturbing acceleration 
resolved along the radial } circumferential, and orbit normal directions. 
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Appendix B 


Retro Motor Optimization Procedure 


The optimization procedure is due to B.A. Glassman et al» 

(Reference 4) and we are indebted to V.C. Zvonkovich of McDonnell Douglas 
Astronautics Company for valuable hints on the implementation of the 
technique. 

The procedure finds values of a given set of independent parameters so as to 
maximize (or minimize) a payoff function while at the same time constrains a 
set of dependent variables. This procedure is based on first-order perturbation 
theory and relies on partial derivatives of the independent parameters with 
respect to the payoff function and the dependent parameters. The partial derivatives 
are determined by numerical techniques. The payoff function is minimized by 
treating it as a constrained dependent variable. The constraint value is systimat- 
ically reduced unt il convergence can no longer be obtained. A derivation of the 
technique follows: 


Let the independent variable be x and the dependent variables be y., then the 
error, e., and the change in the independent variable can be obtained from 


r+1 

x 



(1) 

( 2 ) 


where the bar denotes the desired values 
r is the iteration number and 


[J] 



i - 1, m dependent variables 
j = 1, n independent variables 
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If n is greater than m, an infinite number of x's satisfy eq. (2). An additional 
ci-iterion is imposed on the problem such that the vector, x, is chosen which 
minimizes the step size defined by 

Y\ 

2 

ds =2 




( 3 ) 


In matrix form eq. (3) becomes 


-l 

ds = ( Ax) T D Ax 
. a 


( 4 ) 


where D is a diagonal matrix with 

a 

2 

x. along the diagonal. 


To find the particular solution which minimizes equation (4), first form the 
augmented function 


ib = (Ax) T Dj 1 Ax+ [ e T + (Ax) T J T ]X 


( 5 ) 


where the elements of \ are Lagrange multipliers. The term in the square 
brackets is the transpose of equation (2). Next equate to zero the partial 
derivatives of ?£) with respect to each Ax,. Repeating the differentiation for all 
Ax. and equating all derivatives equal to zero results in 


-1 t 

2D Ax + J X =0 

a 


( 6 ) 


Solving equation (6) for Ax yields 


a 1 D t T % 

Ax = a J x 


(?) 
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Differentiating equation (6) again for Ax. produces 


2 , 

5 ' ' 

3(Ax) 2 



( 8 ) 


Since D -1 is positive definite, the Ax. obtained by equation (7) is a minimum 

a 1 * 

value for i p . 

Substituting equation (7) in equation (2) and solving for \ yields 

X = 2 ( J D J T ) e ( 9 ) 

a 

By substituting equation (9) into equation (7) , the final expression for Ax is 
obtained: 

Ax = -D J T (J D J T ) -1 e ( 10 ) 

aft 

Successive iterations using equation (2) and equation (10) are performed until the 
constraints are satisfied to within an input tolerance. The payoff function is then 
inspected and incremented. 


For the particular application of calculating the retro attitude, the independent 
parameters are the right ascension, 6, and declination, 5, of retro firing attitude. 



where V,^ 


is the trim velocity (payoff function) 


Equation (10) becomes: 



e 2 a.v T 
a e 

5 2 SV T 

a 6 


AV 


T 


2 2 2 2 

e av T + 6 a v T 

3 0~ s S 
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( 11 ) 



Where AV T 


is the error in the payoff function. 


The following strategy is used to minimize the payoff function, 

1. The initial payoff function is calculated and incremented 
by an input amount ( -,ol for RAE-B) 

2. The partial derivatives are determined numerically and 
the equations (11) used to determine the increment in 
right ascension and declination. 

3. The value of the payoff function is determined at the new 
attitude. 

4. The gain in the payoff function is decreased if the payoff 
function is increasing or cannot converge to the desired value. 
The gain is increased if convergence is achieved. 

5. Iterations are continued until the gain in the payoff function is 
within some input stopping criterion . 


Figure 3.24 presents an example of a practical situation. This figure shows 
contours of constant trim fuel required to achieve the desired orbit. The approach 
trajectory shown is a realistic one in that it resulted from a 1 meter per second execution 
error in the first midcourse maneuver. The elements of the lunar approach 
hyperbola are as follows; 


a = -7880.09 km 
e = 1.392407 
cc =139.715 
i = 120.335 
Q = -74.131 


Relative to true lunar equator and 
prime meridian of date 


h m 

Time of closest approach June 15, 1973 5 15 GMT. The path of the steepest 
decent procedure is shown on the figure. The initial attitude used was the velocity 
vector at closest approach. After five iterations, the solution converged to the 
minimum value of trim fuel. 146 



Appendix C 


Targeting Procedure for Inclination 
and Closest Approach 


The impact parameters, B°T and B-R, are components of the vector distance from 
a central body to the point of closest approach of a hyperbolic asymptote for a hyper- 
bolic trajectory relative to that body. They are of value as targeting parameters for 
lunar and interplanetary missions. The more direct objectives of inclination and radius 
of closest approach may be used in defining desired values of B*T and B*R if the 
hyperbolic excess velocity vector is also specified. 

The coordinate reference 'for the impact parameters is 1) the direction of the in- 
coming asymptote or hyperbolic excess velocity vector, S, of the hyperbola, 2) the 
outgoing node of the hyperbola, T, on the equator of the target body and 3) the. orbit 
normal, R. The T-vector is defined by 

T = K x S (normalized to unity) (1) 

where K is the polar vector at the planet. R is then defined by 

R =SxT. (normalized to unity) (2) 

The hyperbolic excess velocity vector, S, is computed from the state vector at 
arrival. Then T and R are computed as Indicated in equations (1) and (2). The 
desired inclination would be achieved (for constant S) by rotating R and T about S 
until the dot product of R with K is equal to the cosine of the desired inclination. 

The sketch shows the S, T, R and K vectors prior to this rotation. 

equatorial plane 


miss-plane 

Sketch: Impact Plane and Miss Vector Geometry 
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Denoting the rotated R -vector by R^, we may write equation (3). 

R = R cos 9 + T sin 9 (3) 

r 

The desired condition is 

K*R = cos i, (4) 

r d 

where i, is the desired inclination. If the declination of S relative to the equatorial 
d 

plane is less than the desired orbital inclination, the combination of equations (3) 

and (4) will provide two solutions for 9 . A choice between these solutions may be 

made based on the sign of i^. This choice affects the node of the resulting orbit. 

The desired B- vector will have no component along R but will lie entirely along the 

- T direction, where 
r 

T = T cos 9 - R sin 9 (5) 

r 

and 9 is the chosen solution of equations (3) and (4). The desired B-vector components 
are thus given by equations (6) and (7). 

(B* T)^ = - b cos 9 (6) 


(B«R) d = bsin9 (7) 

where b is the magnitude of B . This magnitude is derived from the desired distance 

of closest approach, r , as follows. The hyperbolic excess speed is defined by 

P 

equation (8) . 

s = s/s-S (8) 

The velocity at closest approach is provided by equation (9) 

v = 

P 




Equation (10) defines the half-angle of the desired hyperbola 


, sr v . • 

a = tan" 1 ( B_E_) (10) 

Finally, b is defined by equation (11) 

b = ( ^2 +r p ) sin « (11) 

When b from equation (11) is used in equations (6) and (7) along with the solution 0 
from equations (3) and (4), impact parameters have been defined from desired 
inclination and closest approach under the assumption of invariant S. In practice, 

S varies somewhat during an iteration, but the described impact parameter definition 
results in convergence to the desired inclination and closest approach regardless of 
the size of initial errors. 
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Appendix D 

Variable Target Inclination Procedure 


This appendix describes the method for computing desired miss vector components, 

B* T and B»R, for an out-of-plane over-burn. The out-of-plane maneuver is assumed 
to be executed impulsively at periapsis of the approach hyperbola. Sketch B-l shows 
the approach hyperbola, the incoming asymptote, S, the miss vector, B, and the 
periapsis velocity, V , 


Sketch B-l: Approach Hyperbola 



In the sketch, a is the half-angle of the hyperbola. The direction of V may be 

3? 

written as equation (1). 

A A, A 

Vp --B cos ci + S sin ct (1) 

The A symbol indicates "normalized to unity. " Both S and a may be computed from 

A A 

the hyperbolic state vector as may also be T and R, the reference miss vector 

A A 

directions. The miss vector B, may be written in terms of T and R. ... 


B = (B-T) T + (B-R) R 

The direction of B may be written in terms of a rotation, 9, about S measured from 

A a 

the negative T direction toward R. 


* ~ A 

B = - cos 9 T + sin 9 R 

A 

The normal to this approach orbit, H, is given by equation (4). 


A A A A A 

H = BxS = cos 9 R + sin 9 T 


( 3 ) 
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( 4 ) 



If the over-burn situation exists, the retro velocity is fired out-of-plane at an 



The angle at which to apply is computed from equations (5) and (6). 

cos A = <v c 2 - Av 2 - v p 2 ) /v p / Av r / 2 


sin A = ± \j 1 - cos 2 A 
The azimuth change, 8* is then computed from 


( 5 ) 

( 6 ) 


cos jS = (Vp + Av r COS A) /v c (7) 

and 

sin 8 = Av sin A/v . ( 8 ) 

^ r c 

A 

The circular orbit’s resultant normal, H , is then written as 

c 

* * » (9) 

H = H cos fi ± V sin 8 
c P 

or 

/» A A 

H = cos 0 (R cos 8 ± T cos a sin£) 

C A A * 

+ sin0(T cos +R cos a sin £) ±S sin a sin 8 . (10) 
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The condition that the circular orbit should have the desired inclination, i , is 

d 


A A 

H • K cos i, 
c d 


(11) 


where K is the polar vector. 


The algorithm which has been implemented is to call subroutine ORIENT with 


H (9 = 0) = R cos 8 i ( T cos a - S sin a) sin 3 (12) 

C 1 

A 

to be rotated about S until equation (11) is satisfied. ORIENT returns two 9-solutions 
which, from equations (2) and (3), specify the desired miss vector components. 


(B-T) d = - b cos 9 
(B' R) = b sin 9 

U 


(13) 
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Appendix E 


Variable Approach Guidance Procedure 


This appendix describes the method for computing "desired" miss-vector components, 
B*T and B-R, for an in-plane overburn to circularize with retro alone at r^. The 
method involves only computation of the magnitude of B, since the direction is governed 
by the desired inclination as described in Appendix A. The magnitude of B can be 
defined as a function of the radius of closest approach, r . 

r 


, -1 v, r \JC~ + 2 u/r . 

b = (r + li/C ) sin < tan he p» 3 p ) 

P 3 4 


( 1 ) 


In equation (1), is energy, v^ e is hyperbolic excess speed and p is the target's 
gravitational constant. The arctangent in (1) is the half-angle between asymptotes 
of the hyperbola. 



Figure C. 1: Variable Approach Guidance Geometry 

In attempting to solve for r (and also in writing the expression for b), we assume that 

P 

the energy, C , will not change much between trials. Under this valid assumption, the 
"desired" b will serve to pull radius of closest approach in to its desired value. The retro 
velocity magnitude changes very little with changes in closest approach distance and is 
also assumed to be constant in the following derivation of r . Let the radius vector to 

Ir 

the point on the hyperbola where r = r d be denoted by R^ and its normalized vector 

Ak 

by R . This point could occur either before or after periapsis. Let T denote a normalized 
vector lying in the orbital plane normal to R in the general direction of motion. These 

A A U 

vectors R and T, are not to be confused with the miss-plane vectors of the same name. 
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The velocity vector at R can be written 

d 


V = v R + v . T 
r t 


( 2 ) 


where v and v^ are radial and tangential components of velocity at r^. The desired 


post-retro velocity is the circular velocity at r^. 


V = v T 
c c 


(3) 


The required velocity impulse from the retro maneuver is AV. 


AV-V -V=(v c -v | )T-y r R 


(4) 


The problem can be reduced to scalar form by squaring AV. 

2 9 2 2 2 
AV • AV = 5v = v ^ + (v - vj = v - 2 v v. +v 
r c r etc 

If we assume that 5v is the actual velocity impulse expected from the retro burn, 

we can solve (51 for r . 

P 


v = C 


+ 


3 r 


(5) 


( 6 ) 


r 

= _E_ 




(7) 


V = \/^- 
C V r H 


(8) 


5v 2 = C + 

3 r d 




c„- 

3 r 
P 


C, 

3 p 


+ 2mv{*“7= [ 

p i2 yjT l 


2 ^ 3 

6v - C. 

3 r 


]} - 


Equation (10) is quadratic in r . Denoting the curly-bracketed term byfi, we 

P 


can write the solution as (11). 

(i? + c s /- M ]/ c 


r = 
P 
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(9) 


( 10 ) 


( 11 ) 



The sign ambiguity in the quadratic solution was eliminated by the fact that r 

P 

cannot be negative. This value of is used in equation (1) to compute b. 
When the miss - vector errors are nulled by targeting, the radius of closest 
approach is r^ as computed in equation (11). If, as well, the retro impulse is 
applied at as in equation (4), the post-retro orbit will be circular at radius 



Appendix F 


Attitude Errors on a Minimum Fuel Midcourse Correction 


Introduction 

It has been observed that for a particular minimum fuel midcourse correction, right 
ascension and declination errors map into the same straight line when resultant . 
radius of closest approach is plotted against resultant inclination. This result will 
be shown to be characteristic of minimum fuel corrections. 

Attitude Errors 

Small attitude errors cause midcourse velocity variations (vectors) which lie approx- 
imately in the plane normal to the nominal midcourse velocity vector. Denoting the 
targeted midcourse correction velocity impulse by AV and a polar (or other non- 

A. 

colinear) unit vector by K, we can describe the general velocity deviation vector, 
fiAV, which is normal to AV. 

A A A 

fiAV = 01 K x AV + j3 AV x (K x AV) (1) 

A 

In equation (1), a and £ are scalar coefficients of the "east" vector, K x AV, 

A A 

and the "north" vector, AV x (K x AV). These coefficients are functions of the 

A 

attitude errors and | K x AV f in our case. 

Mapping 

Small midcourse velocity variations map linearly into variations in end conditions. 

It is this fact that allows targeting of the maneuver in the first place. Suppose that 
we denote the radius of closest approach (or B*T) by ^ and inclination (or B*R) by 
^ 2 * Let us also denote their sensitivities to AV by 


A = ( /b AV) T 

(2) 

b /a av) t 

(3) 
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In equations (2) and (3), A and B are assumed to be column vectors. Midcourse 
velocity deviations are mapped into end condition variations by equation (4). 



6 AV 


( 4 ) 


We now perform some vector algebra, the reason for which will soon be apparent. 

= A T 6AV 

rp A rp A A 

= aA (K x AV) + £ A AVx(KxAV) 

= - O' K T a X AV + £ (A x AV) T (K X AV) 

( - ff K + flK x^) T A x AV = D T A x AV (5) 

Similarly 6^ 0 = D T B x AV (6) 

A 


Minimum Fuel Guidance 

In what has been done so far, the minimum -fuel aspect has no bearing. It soon will. 
The MFG velocity correction impulse is computed from 

AV = -(7^)" 1 A^ ( 7 ) 

where A 'I' is the original constraint error vector and where (7^) 1 is the 
pseudo-inverse of v ^ . 


v# = 




The pseudo-inverse is obtained here by adjoining a third row, C , ,to 7 ^f. 
C is selected to be the normalized cross-product of A and B. 


( 8 ) 


C = 


= (A x B) / 1 A x B | 


A x B 

The inverse of the augmented gradient is then (by Cramer's Rule) 


( V>£) 




B x C 


C x A 


[ A x Bf f A x Bj 


( 9 ) 


( 10 ) 
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from whence A V = 


B x C 


( 11 ) 


A x Bf 


C x A 1 . , 
A>Jr . 

|Axb|J 


We are interested now in the terms A x A V and B x AV from equations (5) 
T T 

and (6). Since A C = B C = 0, we have 


A x AV = - 


and B x A V = - 


-C (A T B) 

| A x B| 

-C (B T B) 
I A x B| 


A*=k lC 

f A x B| J 

^LU=k 2 c 

| A x B| J 


( 12 ) 


( 13 ) 


We note (perhaps with some surprise) that both A x A V and B x AV are proportional 
to C. If we denote the bracketed vector from (5) and (6) by D as shown there, we may 
re-write those two equations as (14) and (15). 


T 

5 01 = k D C 

6 0 2 = k 2 D T C 


(14) 

(15) 


The relationship we originally intended to prove is (16). 
k l 

6 *l = ~ 6^ (16) 

2 

The attitude variation is represented by the parameters, Q-and £ , of equation (1). 
These parameters then appear only in D, which cancels out in forming the quotient 
of the two end condition variations. 

Geometrical Interpretation 

It may be observed from equation (11) that the minimum fuel AV lies in the plane 
formed by the constraint sensitivity vectors, A and B. 
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C x B 

K 

\ 

\ 



Figure A, 1. Sensitivity vectors in the miss-plane 

As may be seen from the figure, C x A and C x B lie in the plane of A and B, since 

C is out of the plane. AV is a linear combination of C x A and C x B and therefore 

lies in the A, B plane. A, plane normal to AV is also seen to be normal to the A, B 

plane. Variation vectors in the plane normal to AV (epitomized by attitude errors) 

may be projected along A or B. These projections represent the constraint errors as 

formulated in equation (4). The projection of the variation vector along A or B is 

identical to the projection obtained by first projecting the variation vector into 

T 

the A, B plane and then along A or B as seen in (17), because A C = 0. 

A T (I - CC T ) 6 AV = A T 6 A V (17) 

The locus of points in the A,B plane corresponding to small attitude variations in 
AV is therefore a straight line along the intersection of the A,B plane and the plane 
normal to A V . 
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